General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



\ 


V 


J 


IBASA-CB- 1697*4) -HE ISOTBOPIC COBDITIOB OF 
ENEBGETIC PABTICLBS EMITTED FB OH i LABGE 
SOLAB FLABE Pb.D. Th«sis CCAlif°E®A4 iMt. 
ot T*chO 202 p EC AlU/EF AO 1 CSCL 03B 


UdcIas 
024U 


o 


:vK 




1 


h 




s • 


a 




The Isotopic Composition of Energetic Particles 
Emitted from a Largs Solar flare 


These by 
John Spalding 


In Partial Fulfillment of the Requirements 
for the Degree of 
Doctor of Philosophy 


California Institute of Teohnology 
Pasadena. California 


1063 


(Submitted November 23. 1062) 


(SRI 82-11) 


Acknowledgements 


I would Uko to thank ay advisor. Edward Stono. for suggesting this project 
and for providing scientific and personal help and guidance. I would also ilka to 
thank Richard Hewaldt for all the help and oritiolam that he has given me 
throughout my graduate career. 

I was fortunate to be associated with the HIST project. HIST was conceived, 
built, tested, and calibrated by a team led by Edward Stone and Rochus Vogt 
and consulting of William Althouae, Alan Cummings. Thomas Garrard, and 
Richard Mewaldt I have greatly benefited from discussions with all of these 
men. Many others at SRL have contributed to the HIST project, Including 
Brownley Gauld, Neil Hickey, and Marty Smith. Rick Cook has contributed pro- 
portional counters used in the calibration of HIST and provided advice about 
data analysis and interpretation. Neil Gehrels and Keith Krombel as fellow 
graduate students acted as sounding boards for ideas and schemes. 

1 would like to thank my parents. Margaret and Donald Spalding, for the 
encouragement they have given me, and I would like to thank my wife, Donna 
Livant, for her love and for being understanding and helpful as I completed this 
thesis. 

This work was supported in part by NASA under contract NASS- 20721 and 
grant NGR 06-008-160. I gratefully acknowledge tuition fellowships from the 
California Institute of Technology. 


Wa report measurements of icotopo abundance ratios for 5*60 MaV/nuo 
nuclei from a large solar flare that occurred on September 23, 1978. The meas- 
urements were made by the Heavy Isotope Spectrometer Telescope (HIST) on 
the ISEE-3 satellite orbiting the Sun near on Earth-Sun Ubratlon point approxi- 
mately one million miles sunward of the Earth. We report finite values for the 
lsotopo abundance ratios *®C/ l, C. 10 N/ 14 N, ia 0/‘®0, M Ne/ B0 Ne, M Mg/* 4 Mg. and 
n Mg/ M Hg, and upper limits for the Isotope abundance ratios 8 He/ 4 He, ,4 C/ lB C, 
17 0/ 18 0, and 81 Ne/ TO Ne. 

We measured element abundances and spectra tc compare the September 
23, 1978 flare with other flares reported In the literature. The flare Is \ typical 
large flare with ‘Tow" Fe/0 abundance (£ 0.1). 

For ,a C/‘*C. lfl N/ I4 N, ,a Q/ ,a 0, “Mg/^Mg, and “Mg/ "Mg, our measured 
isotope abundance ratios agree with the solar system abundance ratios of Cam- 
eron (1981). For neon we measure M Ne/ w Ne = 0.109 + 0.028 - 0.019, a value 
that Is different with confidence 97.5% from the abundance measured In the 
solar wind by Gelss st a l. (1972) of “Ne/^Ne = 0.073 ± 0.001. Our measure- 
ment for ^Ne/^Ne agrees with the Isotopic composition of the meteoritic com- 
ponent neon-A. 

Separate arguments appear to rule out simple mass fractionation In the 
solar wind and In our solar energetic particle measurements as the cause of the 
discrepancy In the comparison of the apparent compositions of these two 
sources of solar material. 
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Chapter 1 - Introduction 


The Isotope abundances of the light elements in the tolar tytUm contain 
Information that la uttful to tavtral flelda of tcltnca. Thtorlta of 
nuelaotytthaalt and tolar system formation try to explain the observed laotope 
abundanoet. Geologlata can uae obaerved laotope abundance patterna aa 
tracers of various geophysical processes on the Earth, planets, and meteorites. 

A important observation, that has emerged after years of study, is that the 
pattern of Isotope abundances is remarkable homogeneous. To a high degree of 
accuracy the relative abundances of different Isotopes, for a given element, are 
the same in rocks analyzed from the moon, exteorites, or from Peoria, Illinois. 
A very good first approximation to the observed isotope abundance patterns 
states that the relative Isotope abundances for a given element are indepen- 
dent of the source of the sample containing the element. 

A more accurate model for observed solar system Isotope abundance pat- 
terns postulates that the solar system was formed out of an Initially iso topically 
homogeneous ball of gas or "soup". Some smell present day isotope abundance 
differences might be expected In this model, as the result of physical and chem- 
ical fractionation processes operating on the solar system’s initial isotopically 
homogeneous ancestral "soup". 

It is only recently, that convincing evidence has been uncovered about 
inhomogeneity In the initial isotope abundance patterns of the 'soup" out of 
which the solar system formed. In oxygen (Clayton, Grossman, and Mayeda 
1973), and in magnesium and heavier elements (Wasserburg, u apanastasslou, 
and Lee 1979), isotope abundance anomalies have been found that differ from 
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terrestrial abundances by up to 5-6%. This may bo evidence for Inhomogeneity 
la tho initial solar system Isotope abundance pattern. In neon, Eberhardt at a l. 
(1979) have found much larger differences from the common mateorltic com- 
ponents. A neon-E component was discovered which has over 5 times the 
**Ne/*°Ne abundance ratio measured In the "planetary" component of meteor- 
ites. 

In view of the Importance of the isotope abundance patterns in the solar 
system, il is Interesting to study the isotopic composition of the Sun. The Sun Is 
the largest reservoir of solar system material, containing over 99% of the total 
mass of the solar system. Yet. when most researchers quote "solar system" iso- 
tope abundances, they quote abundances obtained from laboratory analysis of 
terrestrial, lunar, or mateorltic material. The Sun Is a qualitatively different 
type of solar system object than the objects (Earth, moon, and meteorites) from 
which solar system isotopic abundance information Is typically obtained. Infor- 
mation about the isotopic composition of the Sun might be used to chock 
theories about the Isotopic fractionation of the solar system’s ancestral ’ soup" 
or to check the hypothesis that the ancestral "soup" was initially lsotoplcally 
homogeneous. However, observational constraints make Isotope abundance 
measurements of the Sun very difficult, and as a result, Isotope abundances on 
the Sun are not very well known. 

TTie determination of solar isotope abundances with analysis of the spectra 
of radiation emitted by the Sun is a difficult task. Electronic transitions in 
different Isotopes of the same element are so close in energy that weak isotope 
lines may be unidentifiable due to overlap from the wings of adjacent stronger 
isotope lines. Rotational and vibrational transitions have been analyzed, but 
require the existence of molecules, and thus a cooler environment than the 
solar photosphere. The atmosphere above sunspots provides the necessary 
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ooolar environment. but thara ara it 111 many lyitamatlo and observational 
problama. In ona analyili, tha ,# C/ **C laotopa abundanca ratio was datarminad 
to an accuracy of 15%, tha s *0/ ,a O abundanca ratio was datarminad to an accu- 
racy of 35%, and tha ,7 0/ ie 0 abundanca ratio was only datarminod to within a 
factor of two (Hall. Noyes, and Ayrat 1978). 

Isotope abundances have bean measured in tha solar wind for noble gases 
(Caiss 1973). but tha measurement of other elements presents formidable 
difficulties. 

Solar energetic particles (SEPs) are a sample of tha solar system material 
accelerated to high energies by tha explosive energy release mechanisms that 
operate in solar flares. The measurement of the isotopic composition of SEPs 
might be useful for tha study of the isotopic composition of tha Sun or to study 
the flare environment that produced the SEPs. The isotopic composition of tha 
SEPs might also provide useful information about nuclear reections that occur 
during solar flares. 

Several technical innovations have recently been applied to a spacecraft 
SEP detection system, constructed from a stack of silicon solid state particle 
detectors. The innovations Include priority systems to select particles with 
charge Zst 3 for detailed analysis, trajectory measuring devices to reduce the 
mass measurement uncertainty for charged particles, stable electronics, 
improved solid state detectors, and heavy ion calibrations. These developments 
have made it possible to measure the isotopic composition of SEPs. In particu- 
lar, the instrument HIST, flown on the satellite ISEE-3 (Althouse at at. 1978), 
can measure the isotopic composition of nuclei with charges 1 < Z < 88 and 
energies 5 MeV/ nuc < E as 800 MeV/ nuc. 

We will report on the isotopic composition of SEPs, observed with HIST, 
from a large solar flare that occurred on September 83, 1978, We measured 
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isotope Abundance ratio* for , *C/* B C, l8 N/ ,4 N, ,i 0/ ,# 0. *Ne/ "Ne, ®Mg/ ,4 Uf, 
and and upper Umita to the abundance ratioe. *He/*He. ,4 C/ '*C, 

IT 0/ **0, and ,, Ne/ B) Ne. We obtain reaulta that are consistent with solar system 
abundances, except in the case of neon, where the solar “Ne/^Ne abundance 
ratio is controversial. The SEP measurement of "Ne/^Ne disagrees with the 
measurement of “Ne/^Ne in the solar wind of Ceiss el a!. (1072). but agrees 
with the meteoritic component neon- A. 

Separate arguments appear to rule out simple mass fractionation in the 
solar wlntf and In our solar energetic particle measurements as the cause of the 
discrepancy in the comparison of the apparent compositions of these two 
sources (solar wind and SEPs) of solar material. 
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Chapter 2 - Hie Experiment 


The Heavy Isotope Spectrometer Telescope experiment, (HIST), wee 
designed to measure the Isotopic composition of nuclei with energies In the 
range from S to 280 UeV/ nucleon and In the charge range with Z * 3 to 2B. A 
general overview of HIST la available In Althouse ef at. (1878), and here we will 
concentrate only on the characteristic! of HIST that are relevant to msasure- 
manta made In a solar dare particle environment. 

HIST la carried on the spacecraft ISEE-3 described by Ogilvle ef at. (1877, 
1878). 1SEE orbits the Sun at approximately the location of a Sun- Earth libra- 
tion point, about one million miles sunward of the Earth. This orbit is an Ideal 
observation point for solar particle measurements because It Is distant enough 
from the Earth that the Earth's magnetosphere does not disturb the solar par- 
ticle measurements. 

The design of HIST can be understood and partitioned as two separate 
tasks, single particle identification and the generation of rate and telemetry 
Information. We will first describe the physical structure of HIST which Is 
relevant to the single particle detection task and give a brief explanation of the 
particle detection method. We will then discuss the generation of rate and 
telemetry information. Finally, we will present a calculation of the Instrument's 
nominal mass resolution. 

Z2. The HOT Particle Detection Telescope 

HIST Includes a stack of eleven silicon solid-state particle detectors 
arranged to form a particle telescope, and the associated signal processing 
electronics. Figure 2.2.1 is a schematic diagram of the telescope, and table 
2.2.1 is a list of the detectors that comprise it. Detectors Ml through TO are 
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•ilicon surface- barrier solid-state detector*. Dat actor* D4 through D9 ara Li- 
drtftad detectors with a central dat action area and an annular gu**d ring 
(•haded in figure 3.2. 1). which la uaed a a an active anti-coincidence ahivid. The 
nominal detection areaa t r the detector* are 508 mm. 9 for Ml and M2. 600 
mm, 1 for 02 and 1X3. and 020 mm.' for detector* !H through D9. 

Single particle identification in HIST ia accompUahed by a multiple parame- 
ter analyaia of the energy depoeltlon hiatory of an ion alowlng down in the sill- 
oon particle detector*. A a la well known (Bertolmi and Coche I960), for a wide 
range of energlea, the charge collected from a reverse-biased silicon diode la 
proportional to the amount of energy lost by a fast ion passing through the 
diode. In HIST the detectors Ml through D8 are connected to charge sensitive 
preampe and amplifier*. The amplifiers' outputs are digitized by a series of 4060 
channel analog to digital converters, with one ADC for each detector. The 
result la a series of pulse heights that are a measure of the energy loss of the 
ion in each detector. 

A particular innovation of this telescope is the use of two position sensitive 
silicon solid-state detectors (Ml and M2), as a hodoscope to measure the incom- 
ing Ion's trajectory. Detectors Ml and M2 are special surface-barrier solid- 
state detectors that have had their metallic contacts deposited as 24 parallel 
strips with a 1 mm spacing, rather than as one large electrode. On each detec- 
tor the metallic strips on one side are perpendicular t j the metallic strips on 
the other side. Each strip has a separate preamplifier and discriminator Thus 
when an ion passes through "matrix detectors" Ml or M2, the position where the 
ion penetrated the detector is known. Figure 2.2.2 shows a picture of one such 
detector. 
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ThMs 2J2.1 - 0 

HOfF Detectors 



Nominal 

Detector 

Norafoal 

nominal 

A0< 

Detector 

Thickness 


Gain 

Threshold 


(pota) 


(lleV/chan) 

(lleV) 

(cm* atr ) 

VI 

60 

sb-Mdet* 

0.121 

0.30 

— 

U2 

60 

■b-Mdet 

0.121 

0.30 

0.79 

D1 

90 

sb k 

0.227 

0.64 

0.73 

D2 

160 

sb 

0.296 

0.71 

0.73 

D3 

600 

•b 

0.617 

1.46 

0.72 

D4 

1700 

U-V 

1.160 

2.76 

0.70 

D6 

3000 

U-D 

1.631 

3.64 

0.66 

DO 

3000 

U-D 

1.531 

3.64 

0.60 

D7 

6000 

U-D 

2.272 

6.41 

0.53 

D6 

6000 

U-D 

2.272 

5.41 

0.48 

DO 

3000 

U-D 

— 

0.19 

0.40 


* surface- barrier detector * “matrix detector" 

* surface-barrier detector 
0 lithium drifted detector 

* estimated error: ± 0.01 
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flfon 2-2.1 - A schematic diagram (to seal*) of the HIST particle detection 
telescope. The shaded areas of detectors D4 to DO are annular guard rings 
used as an active anti-coincidence shield. 
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Agars 22J2 - A photograph of a matrix detector showing the 24 separate 
charge collection electrodes. The electrodes have a nominal spacing of 
« 1 mm. The other side of the matrix detector also has 24 charge collec- 
tion electrodes, which run in the orthogonal direction to the electrodes 
shown. 
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SL3. Single Particle Identification 

1 rtlcle identification In HIST is baaed upon a two parameter analysis 
method called dE/dX - E, that la valid for Iona that atop In the particle detector 
a tack. Figure 2.3.1 ahowa an Ion of charge Z, maaa M. and energy E. which 
entered the teleacope a tack at an angle t), and then atopped In the detector 
stack. The ion haa passed through a detector of thickness T, depositing there 
an energy A£\ and stopped In a second detector, losing there an energy E. If 
the range in allicon, of an Ion of oharge Z and maaa M. Is R(Z,M,E), then the 
equations 

^3- = R(ZM.E) - R{ZM.E) , (2.3.1) 

and E • El 4 LE . (2.3.2) 

express the relations between the variables. 

Equations 2.3.1 and 2.3.2 are two equations for the three unknowns E, M, 
and Z. For a given charge estimate a*, there la a solution M( •*) of 2.3 .1, for 
each Aff and E. Thus, it seems that we must estimate the charge Indepen- 
dently, If we are to measure the Ion's mass M. Actually, there Is a third implicit 
equation that allows a solution for both Z and M. We can use the fact that the 
incident nuclei are stable as a constraint on the number of Isotopes that can 
exist for a given charge Z. For example, a U N ion is Incident on the telescope, 
and deposits 94.62 MeV in D2 and 41.22 MeV in Dl. The solutions are then M = 
20.844 amu for * 8, U = 14.003 amu for a* = 7, and M = 9.921 a mu for a* s B. 
The s* = 0 and * * = B solutions can be eliminated, because the solutions, 
i#(s*), are not near the masses of stable or long-lived radioactive isotopes. 
Then the only solution for a stable isotope is M ~ 14.003 amu, and Z = 7. 
Throughout the entire operating range of HIST, unique solutions for M and Z 
are always possible. 


Figure 2.3.2 is a plot of vs. E for ions that stop in D2. for all stable 
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isotopes with chargea In the Interval 2 * Z « 14 . Bach labeled line, or track, is 
the locus of all possible Ef pairs for a given ion that stops in D2. Here, it is 
easy to see that there is no ambiguity in the solution for 11 and Z. The Isotope 
tracks shown are well separated. 

Figure 2.3.3 is a plot of HE vs. & for data spanning the time interval from 
1078: 286 to 1078:272, for ions that stop in detector D2. The figure 2.3.3 data 
set consists of "raw*' data. No consistency requirements have been placed on 
the data to reduce the background apparent below the carbon track, for exam- 
ple. For the following two reasons, only charge tracks are visible. First, ions 
incident at different angles have different track positions, and the figure 2.3.3 
data set accepted all incident angles out to « 25* . Secondly, each charge has 
only one abundant isotope We will later find, for example, that ,0 C has an 
abundance only IX as large as l *C, so that separate Isotope tracks may be too 
faint to be seen on this type of plot. 
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flgurs 2.3.1 - A schematic diagram illustrating the basic components of the 
dE/dX - E technique of particle detection. 

» 
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flgure 2.32 - Calculated tracks of AE » a. for all stable Isotopes with 
3 ss 2f 14, that stop In detector 1X2. AE Is the energy deposited In the 90/a 
thick detector Dl. E> Is the energy deposited In the 150 p thick detector 
D3. 
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flgora 2.3.3 - A plot of A£ va. £ with experimental data for partlclea that 
■top In detector D2 Each point correeponda to ona event The figure 2 3 3 
data tat conatata of "raw'* data. No conalatancy raqulramanta hare baan 
placad on tha data. LE la tha energy deponlted In the 00 p, t’uck detector 
Dl. A' la tha energy dapoaitad In tha 150p thick detector DQ. 
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8.4. Bmt and Rata Records 

An event rvfsr* to the simultaneous occur ranee of a signal above the 
detector threshold* In detectors M 1 and 1(8. The telemetry record kept by HIST 
of the set of pulse heights hodoecope discriminator triggers and and other 
Information, coincident with an event trigger, la called an event record. The 
deepeet detector In the stack triggered by the particle determines the RANGE 
An event Is a RANGE 0 event If M2 and Ml are triggered An event la a RANGE 3 
event If D3 la the deepest detector triggered. In an event record up to five 
detector pulse heights are recorded. The pulse height for Ml. PHAMl. Is always 
recorded. The pulse height from the deepest triggered detector In the stack, 
PHA1, la also always recorded. Up to three pulse heights from the next three 
detector j above the deepest triggered deter tor are recorded The three pulse 
heights are named PHA2, PHA3, and PHA4, with the convention that PHA2 is the 
pulse height of the next detector above the deepest triggered detector in the 
stack. A particle that enters the stack and stops In D6, has PHA1 for the DC 
pulse height, and PHA2 for the D4 pulse height. 

The pattern of triggered strips In the matrix detector Is not store d directly 
In the event record, but is processed to extract the significant information. 
Any set of as many as four adjacent triggered strips is called a group The 
hodoscope electronics counts the number of strips In a group, and the number 
of groups on each side (plane) of each matrix detector. If only one group of 
strips Is triggered In a matrix detector plane, the event record contains the 
address of the highest numbered strip, and the number of strips in the group, 
called the strip count. For two groups of strips, the event record contains the 
addresses of highest numbered strip and the strip count for each group. With 
two or more groups of strips In any matrix detector plane, we have a condition 
known as the multiple hodo condition. or MH. The MH condition prevents us 



- 18- 


from assigning a unique one perl tele trajectory to the pattern of triggered 
■tripe in the hodoecope. For three or more groups of stripe, the erent record 
contains only the strip address and strip count for the first group, as well as a 
record that three or more groups have fired. 

HIST also classifies an event as one of either PEN. 1DZ or LOZ. The PEN 
class contains all events that trigger detector DO. the last detector In the 
detector stack. All "stopping events", that Is, events that trigger detectors M2 
through D0 as the deepest detector, ire either HJZ or LOZ The HIZ condition is 
a requirement that the charge of the particle he Z a 3. as determined by two 
parameter pulse height analysis. The LOZ condition Is a requirement that the 
charge of the particle be Z <2. This Is acre tip 11 shed with a "box discriminator". 
The HIZ, LOZ, and PEN conditions can be written as a set of logical equations, 
with 


Z3 * ( PHA 1 > Hj) OR ( PHA2 > A^-,) OR (2 4 1) 

[( PHA\ > Lf) AND ( PHA2 > 

HIZ * Z3 fPO 'J/ 1 *U2 (2.4.2) 

LOZ * 73 *U 1 *H2 (2.4.3) 

PEN * D*H\'U2 (2 4.4) 


Here, J specifies the detector, and Hj and L) are the commendable high and low 
discriminator settings, listed In table 2.4.1. 
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TftbU 8.4.1 - 


Ml 

M2 

01 


(ohiDMla) 

(M*V) LO 

(ehuiuii) 


H DtscrtmtaaLor 

CM-V) 

144 

12.92 

120 

10.14 

199 

11 59 

104 

909 

129 

19 49 

99 

959 

139 

23 49 

99 

12.37 

99 

41.33 

59 

1091 

129 

92.94 

104 

59.93 

129 

107.99 

99 

92 33 

129 

109.02 

99 

93.79 

129 

159.97 

99 

90 45 

129 

190.91 

99 

94.14 
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H1ST has 169 rate registers, that art read out to l SEE- 3 every 64 aaconda, 
to maka a rata record Tha rata registers count tha numbar of occurrancaa of 
various coincidence requirements, and ara Uatad in tabla 2.4.2. Rataa 1-16 
rafar to "Motored" rataa. 1SEE-3 la spinning about a fixed axis in apaca, normal 
to tha acliptic plana, and tha aactora a pacify tha 46* octant, ralattva to tha Sun. 
Into which HIST potnta. Rataa 1-6 accumulate only whan HIST la not buay pro- 
ceasing trl 4 f ara or generating avant records, and thus measure tha Uva time In 
aach sactor. Rata raglatara 2V through 169 are only accumulated part time, 
they ara active from 1/3 to 1/24 of tha time depending on tha specific rata 
register. 

Ifcbte 2.4.2 - 

HDT Rata Raglatara 

1-6 60 khz. sectored live tiros clock 

9-16 sectored HIZ events 
17-26 range distribution of HIZ events 

26 LOZ events 

27 PEN events 

28 ADC triggers 

29-124 "natrtx detector single strip triggers 
126- 1 .2 i ange distribution of LOZ events 
133-142 single detector triggers 
143 146 matrix detector triggers by plana 
147-149 single detector triggers 

150-169 guard ring triggers 


I 
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&B. Talatnetry and Event SanpUng 

A very Important characteristic of the 1 SEE-3 spacecraft system la the 
telemetry bandwidth allocated to HIST, 16 eight bit bytes per second, for com- 
munications from the spacecraft to the ground. HIST has the capability to 
court! particles at rates of over 5 x 10 4 particles per second. HIST can also send 
complete event records to the spacecraft at over 6 x 10* particle events per 
second. At the ISEE-3 telemetry rate allocated to HIST, a maximum of only 
* 0.0 event record per second can be returned to Earth. Thus It was necessary 
to design a data sampling system that made some rough cuts In the data, and 
sent only significant samples back to the ground. We will describe a very simple 
functional model for HIST, a. ' then show how it meets the special needs of 
measurements In the solar flare environment. 

The solar flare particle environment In which HIST must make measure- 
ments Is characterized by particle spectra that decline rapidly with energy. If 
F t is the flux per unit energy, for a given element averaged over a solar flare, 
and F t Is fit to the functional form F t ~ AE~^, then y is often found to be In the 
range 3 to 5 The spectral shapes for different elements are often roughly simi- 
lar and the relative abundances are not too different from solar system abun- 
dance tables, such as Cameron's (1961), though systematic trends do appear In 
the data (Cook 1961). 

Suppose there Is a solar flare with /J = \ E~ r *. and all the elements have 
the same spectral form Suppose that the element abundance distribution is 
the same as the Cameron (1961) abundance table. Then, if HIST observes 10 4 
carbon nuclei In RANGE 0, there will be 3030 carbon nuclei in RANGE 1, 464 in 
RANGE 3. 174 In RANGE 3 and 30 carbon nuclei In RANGE 4. If 10 s protons are 
observed in RANGES 0 through 3, then there will be in RANGEs 0 through 3, 
72x10* helium nuclei. 5600 carbon, 6060 oxygen, 616 neon, and 119 Iron 


nuclei. The volar flarv particle environment Is weighted strongly towards the 
low RANGES and hydrogen and helium nuclei. 

What we desire as experimenters, is an event sample with a more equal dis- 
tribution of elements and RANGES. The higher RANGES have better mass resolu- 
tion. and more pulse height measurements of each particle are available as 
checks on measurement errors. When measuring spectra, for a given number of 
events, the least uncertainty in the kpectral shape occurs when there are 
roughly equal numbers of event samples per RANGE. For the measurement of 
the ratio of the abundances of two elements, the minimum uncertainty in the 
abundance ratio is obtained when there are equal numbers of event samples 
per element. 

HIST has an event sampling system that partially meets some of the needs 
expressed in the previous paragraph. HIST has eleven registers in which to 
store event records while waiting for the 1SEE-3 spacecraft to periodically 
request data to send to Earth. Nine of the registers are reserved for H1Z event 
records, with RANGES 0 through 8. One register each is reserved for LOZ and 
PEN event records. HIZ events are given absolute priority. Each time the 
spacecraft requests an event record, a HIZ evsnt record is read out, if there are 
any event records stored in the HIZ registers. If R is the RANGE of the last HIZ 
event record read out then registers for RANGLs R-l, R-2, ... 0, 8, ... and R are 
checked in turn. This feature makes the HIZ event sample have a more equal 
RANGE distribution. Only if no HIZ event records are available, is either a LOZ 
or a PEN event record read out to the spacecraft. 

The measurement process is started when the signal in any one of the 
detectors exceeds the threshold energy listed in table 2.2.1. In a relatively 
quick period of time ( » 10/ii ) HIST decides whether or not to generate an 
event record. HIST first determines if an event has occurred and whether the 
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•▼•at meets th« specified coincidence requirements. Than HIST chacka to aaa If 
there la an empty register In which to ator« th« •want record In tha aaroa 
period of tima tha appropriate rata registers are Incremented. If there la an 
empty a want reftatar of the proper type. HIST win complete the maaaurt mant 
taak. fenerata an event record, and load It Into an event register. This takas 
about 1 me. If thara Is no empty event register of the proper type. HIST wdl 
reset Itself and wait for the next event trigger. Because Hi ST has only tc pro- 
vide event records to ISEE-3 at the maximum rate of • O.w per second it can 
operate at triggering rates *8* 10 4 per second, rather than at maximum 
triggering rates of a 5 x 10" per second, that would be necessary if HIST bad to 
read out every trigger as an event record. 


2d. Theoretical Maas Resolution to HIST 

The solution of equation 2.3. 1 for the mass estimator 

if(s*) * U(9*.L.E.^E) (2.3.1) 

depends on three experimentally determined quantities. E, A E and the path- 

f 

length L * 10 t-h* Aff detector. We will examine how experimental errors 

and physical fluctuations In L, & and Aff affect the width of the distribution of 
the mass estimator. We will follow and expand upon the treatment of Ston » and 
Vogt (1972). 

We will now examine the structure of equation 2.3.1 under the Influer ce of 
two simplifying assumptions. The first assumption Is that the Incident particle 
la going fast enough when it passes through the slab of thickness L that tie Ion 
Is fully stripped of electrons. Then the range of an Ion of mass If and charge Z 

is 


[l. 

wtu,re Ftp Is the range of a proton. Up 


iH* (28z) 

Is the proton's mass, and A = ■— -. A 
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further assumption Is that 




(8.6.3) 


with k * 1 1.71 aieroos of slllooa and X * 1.775. Within ths energy Interval 
r 

10 JfcK* j s 800 U*V aquation 2.6.8 reproduces the tabulated proton range In 
silicon given by Jannl (1066) with a maximum error of loss than 3%. 


Equation 2.3.1 then becomes 




(26.4) 


[*[*-<*)*]] 

The derivatives that we need for our error analysis are given below, 


(2.6.5) 


ul\ *L=lJ* 

OL \g |x-l J L 

OM I . oy_\ _ f x ) u \ r_ 

0g\g 9AE Ix-lJ^lL" 


( 8 . 6 . 6 ) 


(26.7) 


( 2 . 6 . 8 ) 


<* ## > 

where R * R(ZJLE). The uncertainty In U from fluctuations In N independent 


quantities, x<. each with uncertainty e^, Is 


where 


M-ii-r 

■-••• eu 


( 2 . 6 . 10 ) 


( 2 . 6 . 11 ) 


If we look at the structure of equations 8. 6. 6*9 they can be written as 
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Sr ■ < sr> f ' < f • <* 4 w) 

where /« la of order 1. Thus the relative precision Lh*t wo will get In the mui 
measurement la about the Mina aa tha relative precision to which L, l! , and b£ 

ara known aa 


•± 1 , 

U % A-l Ji 


(2.6.13) 


Thua a maaa uncertainty of 0.20 amu at •Ne. aa will bo do mo nat rated, me ana 
that all quantities muat bo known to a precision of bettor than IX. 


2jM. Pathleogth Variations 

Equations 2.6.5 and 2.6.10 show fiat 



1 

X-l 


o_L_ 

L 




(26.1.1) 


One caueo of the variation In the pathlength Is variation in the thickness of the 
LS detector from place to place on tha detector. Since the position at which a 
given Inclilent particle hits a particular detector can be extrapolated with the 
position Information obtained In the Instrument hodoscope, the Instrument 
response can be corrected for thickness variations in the detectors by the use 
of thickness maps. 


The thickness variations of the HIST detectors were measured using parti- 
cle beams. A beam of constant energy was incident on a particle hodoscope 
and then penetrated the detector to be mapped. The thickness of the detector 
can then be calculated using the amount of energy deposited in the detector 
and the value of dE/dX, the energy loss per unit distance, appropriate to the 
particle. The particle hodoscope was an argon proportional counter. Three 
different particle beams were used, protons from the Caltech Kellogg Labora- 
tory van de Graaff accelerator at 8-12 lieV, argon Ions at * 900 MeV/nuc, and 
iron Ions at * 600 MeV/nuc from the Lawrence Berkeley Laboratory Bevalac. 
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ST 

Most of the detectors hod measurable thickness variations, with of ordor 

T 

IX. Tho t hie knots maps won atond on 1 moo. grids, to a precision of * 0.2X. 

A pooolblo limitation of this mapping technique Is its inability to discrlm- 
inato botwoon variations in detector thickness and positional variations in tho 
charge collection efficiency of the detector. Another limitation of the particle 
mapping technique is that the absolute thickness of the detector can only be 
determined to the same accuracy to which the relationship between range and 
energy is known. 

The other source of pathlength variations arises from variations in the 

T 

incident particle’s angle d. expressed in the equation L * ■ „ . Suppose that 

0019 

a particle intersects detector 1(1 at (x,. y ,) and M2 at (x a , y t ), as illustrated in 
figure 2. 6.1.1. For ease in calculation, let y t * y a . and assume that x a is known 
without error. Then. As < x a - x, « Ax + ui, for u> equal to the matrix detector 
strip width, if we assume that only one Ml matrix detector strip has fired. 

*» * w *»•— < 2aiJ! > 

to first order and 


* 7RT (2.e.l.3) 

If both *i and x a are confined to an Interval of width to, on matrix detectors Ml 
and M2, then oi la VS larger. Than we get a Anal result for the pathlength 
uncertainty. 


°L_ 

L 


stnd coed 


w 

l 


(2.6. 1.4) 


with 1 equal to the separation of detectors Ml and M2. This result also holds for 
Vi * y t . At the maximum incident angle for HIST of * 25*. and with parameters 


to 3 i 


Or _ 

and 1 = 50 mm., — = 0.0031. For "Ne this translates to a maximum 
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mw uncertainty ag * 0.001 a mu, For a better estimate we can average In qua- 
ff i 

drature -j- ower an Isotropic distribution of Incident particles so that 


WL_ 


f coed dA dO 
f coed dA dQ 


(2.0.15) 


The case with the largest < —> would be range 0 particles. If we approximate 


Ml and M2 as disks of radius 1.2 cm, and integrate numerically, 

corresponding to a M ( *Ne) « 0.046 amu, and 


<~> * 0.00174, 


on { **Fe) » 0.120 a mu. 
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Wl MO an Ion paw through tha AJT da t actor Lha amount of anarg y loat. AT. 
Auotuataa essentially bacauaa tha nurobar of alactrona that tha Ion "hits'* Is a 
statistical procass. Rossi (1002) dose ribas this energy loss fluctuation process 
from a aemi-clawtcal viewpoint. for thin da tact ora. For thick dot actor* a 
correction factor for tha thin detector result, developed In Spalding (lt'fll). Is 
needed. Tha correction factor accounts for tha fact that a particle's <inargy 
changes significantly as It passes through a thick detector. Tha result Is 

oU * o! k{ff) L If (2 6.2.1) 


o/ = 2* ( 0.136 Jit V )* / (mm. of wiUcon) 


2 ) 


and 


• (S:.6 2.3) 

If Is tha correction factor for the thin detector result and In our poear law 
range energy approximation. 

***** 

D 

for large — , we have a thin detector pathlength and If « 1. Table 2 8.2.1 shows 

ff 

D as a function of — for 1*7/4 and demonstrates that the correc.ion Is 

Li 

Important for y<2. 

u 

IbMe 2.82.1 - 


Thick detector correction factor 


R/L 

D 

R/L 

D 

1.01 

6.33 

200 

1.19 

1.00 

2.77 

3.00 

1.10 

1.10 

2.14 

5.00 

1.06 

1.20 

1.59 

10.0 

1.02 

1.60 

1.34 

20.0 

1.009 
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To eompute tho mui estimator uncertainty wo nood tho derivative of M 
with roopoot to Aff, equation 2.6.7, at fixed E. Thie la because tho energy lose 
fluctuations of A B and E are corrolatod, with 6 (US) * -d (E ), as long as tho 
particle still stops within tho E detector. Then 

•a u • fa f (j - i]* a. (a(A) ^ D <«••*■»> 

and 


»».« * f [(f)^ -(f-D^d-^] • <*•♦»> 

ft 

Figure 2.0.4.1 1 shows ojr.as aa a function of -j- for a 3 mm. detector. As a func- 

R R 

tlon of -j-, Os, es has a finite Intercept, Of.jf(O). at — * 1- and then monotoni- 


cally Increases. For the powor law approximation tho Intercept Is given by 

11 


A 

M 


A— 1 

E 

3A-2 


(2.6.26) 


and in silicon for A * 7/ 4, 

°M.am(0) * 1-712 Z U ( 1 
where E(L) Is the energy for the Ion to go a distance L. For a 3 mm., detector 
this works out to ffjr. as(0) *0.117 &mu for M Fe, and e#. jf(0) * 0.042 amu for 


II 


MAX 


lmm, 


(2.8.37) 


"No. Thtse are fundamental limits for this thickness of detector. 


M.3. Ene r gy Measurement Errors 

The effect of uncertainties in the measurement of Af or A' can be com- 
puted using equations 2.6 8-9 and 2.6.10. In particular, note that at R/L = 1 
errors In A E and if are equally Important since 

p&L'A'ifr <283i) 

where E(L) is the energy for the Ion to go a distance L. 
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Digitisation of the energy signal oauaas an uncertainty where o is tha 

anarg y width of ona paisa haight channal. Except for very snail signals such as 
for proton and helium Ions, this contribution to tha mass uneartalnty is nagllgl- 
bla. For a *°Na ion in range 2 with R/L * 1. c* (dlgiMsoMon) * 0.018 imu. 

An alactronle calibration of HIST was performed, introducing a signal of 
known ebargo at the prtamp inputs by allowing a tost pulM of known height to 
charge up a teat capacitor. Tha calibration was dona at several temperatures 
over tha full dynamic range of the ADC's. Thirty points par detector wore taken, 
approximately logarithmically spaced in energy. Tlie result is that tha charge 
to pulse haight conversion is known to batter than ona half of a channal. Com- 
parison with tha digitisation error shows that electronic calibration errors are 
of the same order and thus negligible. 


A single power law representation of the range energy relationship with 
A/ 2* scaling is only an approximation for the range of energies and ion species 
measured by HITT. To calculate the mass uncertainty, two approximations were 
made, that the proton range energy equation could be approximated by a power 
law form and that ions were completely stripped of their atomic electrons while 
slowing down. At the lowest energies in the HIST operating region both of these 
factors become increasingly Important. We have therefore carried out calcula- 
tions to examine to what extent the mass resolution calculations depended on 
the above two assumptions. 

To evaluate 

o M .u * (fc<0 l]“ D (2.6. 4.1) 

m jkA 

for arbitrary range energy relations, D and zrrj; must be computed numeri- 
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oally. 777 via computed by finite differences end D wee calculated using 
oOE 

numerical techniques from Spalding (1901). 

To aaeeaa the effects of the power law approximation, the range energy 
equation for protons was obtained by Integrating (Marmier and Sheldon 1969) 

dj? ■ | " i | 2" [n * W /° ) * ) " Nl”^ ) “ 0* J . (8.0 4.2) 

with 1 * 170 eV. In figure 2.0.4.1, we plotted Ok.am R/L for "*70 Ions, as a set 
of solid lines, mmunalng complete stripping of the Ion. The dashed lines are the 
single power law approximation of equation 2.0.2 0, and the lines are Labeled by 
the pathlength L. This power law approximation Is a eery good approximation, 
except for R/L a 4 for the 3000^4 detector. This corresponds to R a 12 nun and 
E a 200 If eV/ nuc, or RANGE 7 and 0 events, which are not Important for solar 
flare particles. 

To describe the effects of incomplete stripping, ions were assumed to lose 
energy at a rate 

£-<*•>■ s§L. <*•♦»> 

where Z* Is the effective charge, and the relation 


rl B 

Z 


where q * 1.3* 10“" 


1 


(2.6.4.*) 


f i 

(Barkas and Berger, 1904) wes used. The range of an Ion was then obtained by 

tLF 

Integrating (37)'*. an. am was then calculated for "*Fe Ions, with pathlengths In 
dLX 


the A E detector of 50 fi, 90/4. 150/4, 500/4, and 3000/4, and plotted In figure 
2.O.4.2. The 3000/4 power Law results are plotted for comparison. With Incom- 
plete stripping, the ton does not slow down as fast as when stripped. Thus an 

effective range energy power law exponent " A** = E * would be lower than 

□Us 

for a stripped Ion. A 'ower "A” will raise <Jm.am< as it Is proportional to 
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Also. a slower ion will have a lowor ”X" than a faster ton. Thus, tha thinnar 
datactnra will have largar om.am at tha aama R/L which tha calculation* illus- 
trate. 

In figure 2.6.4 3, a*, a* as. R/L is shown for a pathlength of 50 /* and Ions 
*Ha, l *C, "Ne, "Si. 40 Ca. and "Fa, along with power law results for comparison. 
Hers we can see that only for ions with charges Z a 14 will tha power law 
approximation to a*, at be a serious underestimate for moderate R/L* 2. For 
the thicker detectors the power law approximation will work better. For Z * 14, 
and detectors with L a 00/*, tha power law approximation provides a result with 
acceptable error. Stoce iron is the only element analysed In this report with 
Z a 14, figures 2 6.4.2-3 and the power law approximation span tha necessary 
range to calculate ag ut for this work. 
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fl|MW U.4.1 • A calcu l ati o n of Of n for **r* Iona, assuming complat# 
Gripping. Tha solid Unas usa tha ranga snargy curva obtainad by intagrat- 
Inf aquation 2.0. 4.2 and tha daahad Unas ara tha powar lav approximation. 
Tha Unaa ara labalad by tha aalua of tha pathlanfth L. 




flgurw L6.4J - Calculation of aa.k* for • t F‘« with tha Barkaa and Bargar 
(1964) affactlva charga. Tha Unaa ira labalad by tha walua of tha path* 
langth L and tha daahad Una la tha powar law approximation. 



(nujo) 3y * w jo 
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flgura 2J.4J - Calculation ol pjj.as f or a 50/i pathlangth L and with Barlcas 
and Bargar affaetiva charga for various ions. Tha daahad Unav art tha 
powar law approximation. 



R/L 
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Just as a reduced ”X" lands to increaee . a raducsd "X” will also 

0 u 

Ineraasa xr- The correction for lnoomplato stripping at low energies affects 
OL 

0i# 

tho derivative — used to calculata tha affact of pathleng th variations on tha 

mass astlmator. In tha single power law approximation tha value of this deriva- 
tive given by aquation 2.6.5 is a constant, as a function of R/L Tho derivative 
was calculated numerically using tha Barkas and Barger (1904) affective 
charge. Tha reduced derivative 



numerical calculation) 
OL 

014 

—-(paiumr law a p prox im ation) 
OL 


(2.0.4.6) 


was calculated for "’Fe ions and pethlengths of 50/a. 90/a, 150/a. and 500/a of sili- 
con. and is shown in figure 2.6 4.4. For the 50/a detector, Incomplete stripping 
causes a^.i to be * 50% larger than the value given by the single power law 
approximation, for R/L < 2.30. For lower charges the effect is not so severe, for 


silicon and a 50 /a pathlength. 



1.50 


only for R/L < 1.40, and for 


neon 



1.50 for R/L < 


1 . 20 . 
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flfnr* 2.8.14 - Calculations of ths rsducad dartvativa —r, as dsflnsd in 

OL 

aquation 2.0.4.S, for **Fs and with ths Barkas and Barg or affacttvs charga. 
Tha linos ara labs lad by ths valua of ths pathlangth L 



paonpaj 


1 £ 

m 
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Por Fe Iona, the inclusion of Information about the charge ststs of ths ions, 
represented by ths Barkaa and Bargtr oorrsotlon. in ths calculations for s*. * 
and sjr.av. substantially changes ths results. Por detectors thinner than lhtyi, 
and for R/ L * 3, ths mass resolution is substantially worse than would be 
predicted if the ions were completely stripped of their atomic electrons. The 
Barkaa and Berger correction is only approximate and was derived from meas- 
urements in emulsions, rather than in silicon. We can use the Barkaa and 
Berger correction as a guide to see where the Incomplete stripping of atomic 
electrons Influences the mass resolution. An accurate evaluation of incomplete 
stripping's effects on o^x °mxm- may require more accurate knowledge of 
the Ion's charge state as a function of velocity, than is represented by the Bar- 
kas and Berger correction. 

g fr if Smxntry 

For table 2 8 5.1 we have calculated the nominal mass resolution of HIST 
for RANGE 1, 2. and 3 events in an idealised situation. The detector used was 
the deepest triggered detector in the stack and the &£ detector used was the 
second deepest triggered detector. We used the Barkaa and Berger effective 
charge to do the computations, which were carried out for R/ L * 2. We 

a i 

assumed a flat detector and that — * 0.00174, as calculated In section 2.8 1. 

We calculated the mass uncertainty due to digitisation and ignored the mass 
uncertainty due to uncertainties in the electronic calibration. The mass reso- 
lution uncertainty due to energy loss fluctuations dominates the mass uncer- 
tainty for this choice of parameters Por the heavier ions, the second most 
important contribution to the mass uncertainty is due to thickness variations. 
Tor the lighter ions, the second most important contribution to the mass 
uncertainty Is due to channel digitization. 
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Ikbla 8.SA.1 • 


RANGE 1 - Nominal Man Revolution (R/L * 2) 



(pathiength 

fiudluaUaai) 

(ohamul 
digit leaf Inn) 

(ohttuul 

dlgtUiaUon) 

(total) 

z 

°a.L 

a M.AM 

°M.U 

°M.g 

o a 

2 

aoio 

0.000 

0 062 

0.025 

0.009 

6 

0.030 

0.128 

0.020 

0.014 

0.133 

18 

aooo 

0.186 

0.020 

0.010 

0.100 

28 

0.211 

0304 

0.010 

0.000 

0.430 



RANGES 

- ihninal Man 

Waanhittnn (R/L 

= 2) 



(pathlength 

(™rp loee 
flwtnaUoni) 

(channel 

digitisation) 

(channel 

d^Uaatloa) 

(total) 

Z 

°a.i 

°M.AM 

AM 

°M.g 

Og 

2 

0.010 

0.073 

0.087 

0.023 

0.102 

0 

0.020 

0.117 

0.038 

0.013 

0.128 

12 

0.060 

0.180 

0.026 

0.000 

0.178 

28 

0.107 

0.200 

0.020 

0.007 

0.336 


RANGE 3 

- Nominal Maas 

Raaotution (R/L 

= 2) 



(pathlength 

variations) 

(energy loaa 
fluetuaUona) 

(channel 

(channel 

dlglUeatlon) 

(total) 

Z 

°M.L 

°M. At 

°a. am 

°M.g 

°a 

2 

0.008 

0.080 

0.084 

0.035 

0.101 

8 

0.020 

0.110 

0.034 

0.020 

0.120 

12 

0.057 

0.155 

0.024 

0.013 

0.167 

26 

0.150 

0.252 

0.019 

0.010 

0.294 
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Chapter 3 - Data Analysis 

9ll. Hm Measurement Overview 

VP 

In ths two parameter - f analysis taohnlqua of sootlon 2.3. a single 

mass astlniator was generated. For ths purpose of making accurate isotope 
msssuramants. throe and four parameter measurement techniques will be used. 
Hth three and four measurements of a particle's energy deposition history in 
the particle telescope, several mass estimators can be calculated and com- 
pared. With this typo of redundancy, we can detect and remove sources of 
background events from the data set Also, to reduce the mass measurement 
uncertainty, the measurements from several mass estimators can be used in a 
weighted mean mass estimator. 

Figure 3.1.1a is a sch imatic diagram illustrating a three parameter RANGE 
3 event. A particle has stopped in detector D3, depositing an energy E! and los- 
ing energies LE X and &£g in the two previous detectors. In this configuration 
there are actually three separate ways to form a two parameter mass estimator. 
The "A E detector" can be taken to be Dl, or D2. Also, D1 and DS can be treated 
together, as one detector. Only two of the mass estimators are Independent; 
therefore, once the total energy A * £ + LE X + Aft* 1* computed, there are 
only two parameters left in this system. 


f 

t 


1 
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figure 3.1.1a - A vchematlc diagram oi a three parameter event. A£j. Alt. 
and g are the energies deposited by a particle slowing down In the Indi- 
cated detectors. 

figure SLl.lb - A schematic diagram of a four parameter event. 



t^x 

A*, 

e ' 
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For this three parameter system. we will compute two man estimator*. 
The first estimator, called U n , uses D2 as the "LE detector" and is found by 
solving 

L « R{M m ,Z. LE , ♦ g ) - R{U m ,Z.£ ) . (3.1.1) 

The second. Hit. vises D1 as the ’ A E detector" and can be found by solving 

L * ff(Jf, t ,Z, Aff§ ♦ Aflj ♦ A* ) — *(lf,, ,Z. A*, ♦ g ) . (3.1.2) 

For the special can in which all of the man uncertainty is caused by energy 

bon fluctuations, and in which the fluctuation in AJ5*, d(Afli), l* small compared 

to AJTf. the man estimators J (fa and if ia are statistically Independent. For the 

four parameter case illustrated in figure 3.1.1b we compute analogously Jf M . 

y m and U m To avoid confusion, detectors Ml and M2 are labeled as A and B. 

respectively. Thus the mass estimator that uses detector Ml as the A E detector 

and M2 as the stopping detector is called Ham For all RANGES except RANGES 0 

and 1. three mass estimators are produced. For RANGE 1 we calculate the two 

man estimators U § j and Afei. and for RANGE 0, we calculate ifj*. 

In the following sections, we will describe the calibration of the mass esti- 
mators and then un mass vs. man correlation plots to describe lata selection 
and background reduction procedures. 

32. Man Estimator Calibration 

In order to measure accurately a particle's mass by the two parameter 
tlE 

— - E technique, the range energy relation used in equation 2.3.1 
aJL 

L - R{U,Z,E) - R(U.Z.E) (2.3.1) 

must be known to sufficient accuracy. If d(A R) is the systematic error in 

A R - R{E) - R{1 0. then the systematic error in the mass estimator is 

dif 

M 

Unfortunately, published tables of the range energy relation?, are insufficiently 


1 g(AR) 
A-l A R 


(3.2.1) 


accurate for our purposes and bar* systematic errors of the order of a few par 
cant. We need to know the location of each Isotope track to a 0 2 mu if track 
location errors are to make only a small contribution to the mass uncertainty. 
For **Ne this means that the systematic error in A R must be less than 0.76X. 

One approach used by Vidor (1076) and Cook (1961) la to attempt to con- 
struct a new corrected range energy relation 

R{M.Z.K) - R,(MJ.f) ♦ RccmiU.Z.D (3.2.2) 

where Rom 1* • function of several unknown parameters } a, J The values of 

the a* are found by minimising the distance, in some sense, between the exper- 
imental and theoretical Isotope tracks. 

The Rpm method needs enough experimental points to locate accurately 
the experimental isotope track in a sufficient number of energy bins in each 
RANGE. If 20 particles are sufficient to locate the isotope track in each of 26 
energy bins per RANGE, then one needs about 500 particles per RANGE to cali- 
brate an element Unfortunately, for the set of HIST flight data, the only ele- 
ments with charge Z > 3. in RANGES 0-4, for which there is su an amount of data 
approaching the above amounts, are carbon and oxygen. 

A factor that limits the accuracy of the Rem function is the possibility of 
small systematic errors in the energy calibration or in the thickness maps. 
Since the energy ranges for all the possible mass estimators .1^ overlap, 
attempting to fit one Ram relation to the data for several mass estimators may 
not be completely successful. 

Oxygen data from RANGES 0 to 3 and mass estimators Ujji, U A \> M it . M a> , 
idg \ , JVjb, Mm- if »s. if .a end M a were used In an attempt to find one unique Ram 
to represent the oxygen data. A 16 point Ram function and the thicknesses of 
detectors Ml, M2. Dl, D2 and D3 were the parameters to be fit to the data The 
resulting Ram function was then used to recompute all the mass estimators. 
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Dsrtations of up to 0.6 amu from the correct ,# 0 mua war* found. For ,a C there 
was a raault. Thus a corractioo aohama was developed to oorract each 

mass astlnmtor, separately. 

In aaoh RANGE tha preliminary mass astlmator was plottad ▼» El, tha 
anar gy In tha stopping detcc»or Figure 3.2.1 showa an example of this type of 
plot, for hail urn In RANGL 0. Tha approximate form of vs. El for the princi- 
pal isotope of an element was then fit by hand to the plot by a series of line seg- 
ments. This approximate form of the preliminary mass astlmator, g(El), was 
then used to correct the mass estimator of each particle with the equation 

M* (corrected) = J ^ (preliminary) . (3.2.3) 

where U, Is the mass In amu of the principal Isotope. 

The advantage of this procedure was that It was quick and easy to imple- 
ment on a small computer system. One disadvantage is that the procedure is 
not functionally exact. If the range energy relation Is a function of E/ll, and 
g/Jl varies slowly with El. then My (corrected ) Is equal to M, plus terms of 
order (p/ A^)*. Another ’’disadvantage” of the procedure might be the practice 
of doing the fits to mass vs. El by hand. A difficulty is that you do not get any 
information about the "goodness of fit” or errors in the fit parameters. In prac- 
tice the hand process Is quicker, much less confused by noise and background 
and gives quite acceptable accuracy. The maximum error in this track location 
procedure is about 0.2 amu. 
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Qgura 3JL1 - n El for RANGE 0 hall urn data. Tha aoUd Una U tha hand 
flt. f(El), uaad to oorraot tha prallmlnary maaa aatlmator, using aquation 
3.2.3. 
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aa. HihhMkittalkUSatariMtkB 

To measure Isotope abundances. a high raaoluUon data sat was prepared. 
Wa will doaorlba selection procaai to r that high raaoluUon data aat Restrle* 
tiona on tho hodoaoopa atata, anargy dapoaltad In Dl. radial poalUon In oartain 
datactora, and masa estimator consistency, will bo described. 


To compute a masa aotlmator, tha pathlongth ot tha parUcla In tha A K 
dataotor must ba calculated. Tha basic raqulramant Is that tha hodoaoopa 
informaUon specify a unique parUcla trajectory through tha telescope. Thus, 
wa must have only ona group of triggered strips par matrix detector plana, or 
tha MH condition must be falsa. Tho MH events ware rejected, as ware event 
records In which over two strips par group ware triggered. Tha fracUon of 
event records with tho MH condition was as large as 1 /4 of tho total number of 
events, at times of high parUcla counting rates. Host of tha MH events were 
thought to ba caused by accidental coincidences at times of high singles rates, 
aa detailed In Spalding (1063). Records with over two strips per group were an 
Insignificant part ( < IX) of the data. 

The remaining event records were divided Into two groups, event records In 
which only one matrix detector strip triggered per hodoscope plane, and event 
records In which one or more planes had two adjacent triggered strips. If hodo- 
scope information were all that was required from detectors Ml and M2, since 
tha particle had stopped deep in the detector stack, then both classes of events 
ware accepted. Because the parUcla had to penetrate the detector right 
between the strips, we think that the double strip groups may have even located 
tha particle better In the hodoscope, thereby reducing the position uncer- 
tainty, although this idea has not been put to a conclusive test. 
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The significance of the double strip condition Um in what it tails us about 
tha measurement of the energy deposited in a matrix detector. Suppose that a 
particle has triggered two adjacent stripe on one plane of a matrix detector. 
Then there Is a large probability that a "signal defect" in that matrix detector 
has occurred. The energy measured by HIST In that matrix detector will be lese 
than the energy actually lost by the particle slowing down In that matrix detec- 
tor. There la also a population of events that have energy measurements 
characteristic of the "signal defect" condition, and haws only single strips trig- 
gered. 

The "signal defect" behavior can be Illustrated with a mass vs. mass corre- 
lation plot. Figure 3.3.1. 1 shows M§ x plotted vs. U^ x for carbon with the single 
strip hodosoope requirement. On this plot, if there was no signal defect 
mechanism about 99% of the data should be in a clump at the point (12,12) and 
1% at the point (13,13). The clump at (12,12) has a horizontal streamer of par- 
tldee with llg\ ** 12 and ((u with some arbitrary lower value. This is caused by 
a signal defect in the matrix detector 1(1, used as LE detector for Jfj t , and not 
considered for U§ x . The slanted streamer falling from (12,12) is caused by a 
signal defect tn the matrix detector M2. M2 Is used as the tiE detector for Ug\. 
and as part of the E! detector, together with detector Dl, for i(i x . The streamer 
than slants because the M2 signal defects introduce correlated errors into M§ x 
and jtfg|. 

In figure 3.3.1.2 and 3.3. 1.3. we plot carbon RANGE 1 events in which dou- 
ble strip events are required in the "y plane" of detectors Ml and M2, respec- 
tively. The double strip requirements have greatly increased the signal defect 
probability. Figures 3.3.1.4a-b are histograms of i(| t and Ifa with a single strip 
hodoscope requirement. Figures 3.3.1.5a-b are histograms of U tx and U§ x with 
a double strip requirement in the "y planes" of matrix detectors Ml and M2. 
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respectlvely The man resolution for the double strip eventa is ao bad that wo 
must dlacard thorn from tho high resolution data aot.Tablo 3 3 11 ahowa tho 
fraction of H1Z evonta In tho time period 260: 10 to 272:1)0 In RANGEa 0 through 4 
for oach of four coincidence equations M2d and Mid refer to the double atrip 
condition in datectora M2 and Ml respectively. In parentheaea la the ratio of 
the number of carbon eventa for tho coincidence equation to the number of 
oarbon eventa In that RANGE. 


Table 3.S.1.1 - 


Carbon Double Strip Eventa 


RANCE 

(U2d) *(U\d) 

(M2d)*(JTld) 

(JT2d)*(J/ld) 

{Xfid) Rifl'd) 

1 

130 

666 

356 

1681 


(.060) 

(.266) 

(.138) 

(608) 

2 

42 

262 

184 

1061 


(.027) 

(.164) 

(120) 

(689) 

3 

21 

179 

146 

868 


(.017) 

(147) 

(120) 

(716) 

4 

1 

28 

24 

199 


(004) 

(HI) 

(.096) 

(.790) 


In appendix A we And that the probability of a double strip event In a 
matrix detector la a function of the energy deposited In that matrix detector 
For the Ions analyzed In this report, the probability of a double strip increases 
with the amount of energy deposited In the matrix detector. This effect can be 
observed in table 3.3 1 1 . The carbon tons that stop In RANGE 1 are going slower 
a a they penetrate detectors Ml or M2 than carbon Ions that stop in RANGES 2-4. 
end thus deposit more energy In detectors Ml or M2. The RANCE 1 events then 
have a higher fraction of double strip events than the RANGE 2-4 events. 



- 50- 


We can obtain some Insight Into tha double atrip events by considering 
briefly charge collection In a "normal" (continuous electrode) solid state detec- 
tor of the same nominal thickness and applied voltage as the matrix detector 
Ml. With an applied voltage of 25 volts, electrons would take about 1 3 nsec, 
and holes about 2.5 nsec, to croes 50 p of silicon. The time constant of the 
preamps Is about 2 /* sec, so that In a normal detector, all the charge would be 
collected. In a matrix detector, the strong association of a large signal loss with 
the two strip events suggests that something happens to the charge collection 
procedure In the gap between the two strips. The electric fleld configuration In 
the gaps might cause some fraction of the eh t. irons and holes to be directed to 
regions of low fleld strength. There the electrons and holes may either recom- 
bine or travel to the strip electrodes so slowly that the signal is reduced, 
because part of it Is effectively filtered out by the preamp time constant. 
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flfura 3.3. 1.1 - U 9l vs. M 4l for carbon data. The hodoacopa triggers are 
restricted to single strip events. 
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flgura 13.12- Ubi vs. U Al for carbon data. Tha hodoscopa triggers ara 
rastncted to events with double strip triggers In tha "y plane'* of Ml. and 
single strip triggers in tha remaining matrix detector planes. 
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fl|un 3.3. 1J - tig | v». U 4l for carbon data. Tha hodoacope trtggera ara 
raairtctad to avanta with doubla atrip trigger* In tha "y plana" of M2. and 
dngle atrip trigger* In tha remaining ms ' -lx datactor planaa. Tha atraamar 
of avanta around lUi * 18-21 and tig i * 8-20 la dua to nitrogen avanta 
which have had their maaa aatlmatora calculated with Z=8 Inataad of Z=7 
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figure 3.3.i.4*-b - Mass histogram* or carbon RANGE 1 avents with single 
strip triggers. Figure 3.3.4a shows U^\ and figure 3 3.4b shows U§\ 
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0|art 3.3.1. 5a-b - Maas histogram* of carbon RANGE 1 events. Events In 
figure 3.3.5a hava doubla atrip hodoacop* triggers In detector Ml's "y 
plana” and single atrip triggers in the remaining detector planes. Events 
in figure 3.3.5b have double strip hodoscope triggers in detector M2's "y 
plane" and single strip triggers in the remaining detector planes. 


s 
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3.3.2. Radius Requirements 

The particle trajectory In the bodoacope la extrapolated to all the detec* 
tore that the particle hita. In two detectors, the radial dlatancea from the 
center of the detectors to the positions where the extrapolated trajectories 
intersect the detectors, were confuted. Dl, the smallest detector exclusive of 
the matrix detectors, a^J the next detector deeper In the stack beyond the 
■topping detector, D(R+l), were chosen. Particles that had an Impact point 
nsa^ the edge of a detector were more likely to show a signal defect than others 
nearer the center of the detector. For the high resolution data set we set max- 
imum radii in the two detectors based on the behavior of the mass estimator 
near the detector edge. 

Figure 3.3.2. 1 shows the weighted mean of the RANGE 1 oxygen mass esti- 
mators plotted vs. radius in Dl. At the largest radii one can see a systematic 
underestimate of the mass caused by the signal defect at the detector edge. 
For Dl, a maximum radius of 13.8 mm was set. The oxygen mass estimator 
weighted sum was plotted vs. radius in D2 in figure 3.3 2 2 Again there is a sys- 
tematic mass underestimate and a maximum radius of 15.0 mm was set here. 
For RANGES 2-4 no maximum radius was set in detector D(R+ 1), For the helium 
data set more restrictive radius restrictions were needed, because the rare iso- 
tope *He lay below the more abundant Isotope 4 He. Correlated mass errors pro- 
duced bv signal defects at the detector edge might then masquerade as 9 He 
events. Table 3.3.2. 1 lists the radial restrictions for helium and elements with 
charge Z > 8 for the high resolution data set. 
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Muimum Radii ( mb ) 
Zk| 


RANGE Radius Dl 

1 13.8 

2 13.8 

3 13.8 

* 13.8 

6 13.8 


RuUus D(R+1) Radius Dl 

15.0 13.9 

non* 13.0 

nona 13 0 

non* 13.0 

non* 13.0 


11*11 inn 

Radius D(R+1) 

1290 

16.75 

18.75 
non* 
non* 
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figure 3J21 - The RANGE 1 oxygen weighted aum mass estimator la plotted 
re the rediui in detector Dl. The maximum D1 radius chosen for the tugh 
resolution data set Is shown as a dashed line. 




ORIGINAL PA -'- I- 1 

OF POOH gUA^Ilf 

- so * 

Ogun 33 2 2 - Th* RANGE 1 0 «yl*n «*l«ht*d aura me** •* Um * ,or '* P 10 ' 1 * 11 
„ U,. r4dlul a detector 02. Tb. meetrmm D2 redid, choeen lor th. high 

resolution data sat t* »hown as a dashad Una. 
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X3A Energy Limits 

In four instances particles war# required to pasa a requirement on Cl. the 
energy deposited in the stopping detector For RANGE 1 carbon eve nee we felt 
that the ,a C track poeltion for £ 1 < 7 U»V waa not known sufficiently well, ao 
RANGE 1 carbon events had to have £*1 a 7J UeK to be included In the high reso- 
lution data set Thus we eliminated about 12% of the RANGE 1 carbon data 

For RANGE 1 neon events, U Bx vs U A , cross plots were prepared for five El 
energy intervals of width 28 IfeV In the lowest energy interval, the mass reso- 
lution was much worse than in the other four energy intervals; so therefore 
neon RANGE 1 events were required to have El a 28 JUOK Thus we eliminated 
about 23% of the RANGE 1 neon data. For nitrogen in RANGE 1. the mass resolu- 
tion for £*1 < 30Jl#sV appeared worse than for £1 a 30J#sK. Nitrogen events in 
RANGE 1 were required to have £1 a 30MsK, thus eliminating about 40% of the 
RANGE 1 data 

In RANGE 4, particles that have completely penetrated detector D4 
encounter two “dead layers” which have a total thickness of approximately 100 
microns of silicon, at the rear of D4 and at the front of DO In these dead 
layers, charge is not collected, thus causing a signal loss that systematically 
lowers the mass estimators A signal loss in the stopping detector would cause 
a correlated mass underestimate in both U& and U $ 4 . that would interfere with 
the *He measurement. A particle of a given mass that stops in the dead layer of 
D4 has the same energy measurements in the E and A£ detectors. D3 and D 4. 
as a lighter particle that stops in D4. but in front of the dead layer. Thus, the 
mass measurement Is ambiguous near the "end of range” position in D4 Only 
tor helium is there a significant population of solar flare events in this RANGE 4 
"end of range" position (* 30 out of 207 range 4 events) We therefore chose a 
requirement £1 < 50 U*V, for RANGE 4 helium events, which removes the region 
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of dead layer mass ambiguity. 


Because of the presence of occasional large errors In the mass estimators, 
which would bias the weighted mass estimator, data included In the high resolu- 
tion data set had to pass consistency requirements on the mass estimators. For 
elements C, N, 0, Ne. and Mg, we chose as a mass consistency estimator the ratio 
of two mass estimators. For the three parameter systems the ratio chosen was 




where J Is the detector number. For RANGE 2. /?, * 


*is 

¥ss 


(3.34 1) 


For a four parameter 


system, R x is the same as above and Rg was defined as 


For RANGE 3. Rg = 


M»a 

Urn 





(3.3.42) 


The distribution of R x and Rg Is approximately Grusslan. For each RANGE 
and element the entire HIST data set was used to form histograms for R x and 
Rg. The histograms were then flt to determine the location and width of the 
peaks In R x and Rg. Figure 3.3.4. 1 shows a histogram of R x for RANGE 1 oxygen 
events. For C, 0, and Mg, mass consistency ratio requirements were set, which 
were symmetric about Rj = 1.0, and approximately Included the central two 
standard deviations, 2 X*. of the ratio peak. For N, the ratio limits were set 
exactly at the fit peak location, plus and minus two standard deviations of the 
ratio peak. With Ne, the same procedure was followed but the requirements 
were relaxed to plus and minus 2.5 X^, because we wished to include as many 
Ne as possible, in order to reduce statistical uncertainties in the tr Ne/ 80 Ne 
abundance ratio. 
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Tbt' fJ.tNGE 3 neon consistency requirements war* widened to aceapt all of 
tha RANGE 3 noon data except for vary larga mast disagreements Figures 
3.3.4 2a-b show RANGE 1 and RANGE 3 naon maaa croas plota with tha con- 
aistancy raquiramanta drawn In. Othar axcaptlona ara tha RANGE 1 C, N, and 0 
consistency raquiramanta, which ara about plua and minua 1 5 X* In RANGE 1 
tha background problama wara mora severe than in tha othar RANGEa and tha 
C, N, and 0 maaaura manta wara for iaotopaa at or laaa than tha IX level of tha 
mora abundant Iaotopaa. Thua we choaa a more restrictive consistency require- 
ment, selecting on tha ratio histogram only tha central Gauastan peak and 
excluding moat of tha tails of tha ratio distribution. Table 3 3.4.1 contains tha 
mass consistency raquiramanta chosen for tha elements C. N, 0, Na and Mg 



ORIGINAL PA; _ |i, 
OF POOR QUALITY 

- 83- 


flgura 33.4.1 - A histogram of the mats consistency •■limator. 
/?, * Mbi/M A \. for RANGE 1 oxygen daU. 
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figure 3.3.4.2b - A plot of Mm ** Mu for noon °‘NGE 3 events. in the tun* 
period 70 266 10 to 70 272 00. The mui coniu ncy requirements ere 
shown as dashed lines. 
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IkNs 3.3 4.1 - 


Has Ratio Requirements 


aeraent RANGE 

Ratio 

limits 


carbon 1 

Ubi / Uax 

10 ± 0 060 

1.4 

2 

M lt / M n 

1.0 ± 0 070 

2.1 

2 

Uxz/Uau 

1 0 ± 0 075 

2.0 

3 

Urn /!/* 

10 ± 0 060 

20 

3 

Mo/ M^ 

1.0 ± 0 065 

20 

4 

M u 

1.0 ± 0 060 

1.9 

4 

Mi */ Ux * 

1.0 ± 0 070 

2.1 


nitrogen 

1 

Mbx / M ai 

1.000 ± 0 052 

1.5 


2 

Mxu/Mn 

1 000 ± 0 053 

20 


2 

Ux*/M At 

1 000 ± 0 067 

20 


3 


0 990 t 0 044 

20 


3 

Mam/ Mga 

0.994 ± 0 055 

2.0 


4 

m »/ m 2a 

1.000 ± 0 047 

2.0 


4 

M*/U lt 

0 995 ± 0 050 

20 


oxygen 1 

Mbx / U A x 

1.0 ± 0 050 

1.4 

2 

Mxt/Mn 

1.0 ± 0 066 

2.1 

2 

m 12 / m A2 

1.0 t 0 070 

2.0 

3 

Mo / M x 3 

1.0 ± 0 060 

2.1 

3 

JVa/tfn 

1.0 ± 0 005 

2 0 

4 

i/s«/ A/g 4 

1.0 ± 0 050 

2 0 

4 

M */ M „ 

1.0 ± 0.050 

2.0 


neon 

i 

Max / M ax 

1 000 ± 0.055 

2 5 


2 

Mit/Mjn 

1.000 ± 0 049 

2.5 


2 

Mo / 

1.007 ± 0 004 

2.5 


3 

Mo/M , 

1.000 ± 0 051 

3 0 


4 

M*/M„ 

1 000 ± 0 002 

25 

magnesium 

2 

Mlg / Mb 2 

1.0 ± 0 04 

1.0 


2 

Mx»/M az 

1.0 ±0 05 

23 


3 

Mo / Mya 

1.0 ± 0 04 

1.9 


3 

Mo / M n 

10 ± 0 05 

2 1 


4 

Ma */ M z 4 

10 ±0 04 

1.7 



. 67- 


ORIGSNAL PAGE IS 
OF POOR QUALITY 


The choice o( a matt consistency requirement baaed on mass ratios implies 
that is proportional to M for dillerent isotopes of the same charge. We would 
expect such a result if most of the contribution to an came from uncertainties 
in the detector thickness, or in the Isotope track location. On the other hand, if 
most of the contribution to an came from other sources, such as energy meas- 
urement uncertainties, or energy loss fluctuations, then an would be more 
nearly constant. For C, N, O, Ne, and Mg, the variation in M is not sufficient to 
test the behavior of an 

For helium, there are sufficient events to check the scaling of aj t experi- 
mentally. The time intervals 296:09 to 297:09 and 307:00 to 306:00 contain two 
3 He rich flares. RANGES 1 and 2 have enough events to measure ojr for both *He 
and 4 He, separately (see figure 3.4.9). Table 3.3 4. C contains an for that time 
period as a function of RANGE aod isotope. For helium, in this case, the data 

are consistent with the hypothesis that an is independent of M. 

Table 3 . 3 . 4.2 - 

Helium Haas Resolution 

Mass <*|f(*He) oj( 4 He) 

U BX 0.118 £ 0.010 0.102 £ 0.020 

U AX 0.238 £ 0.037 0.228 £ 0.029 

U xi 0.138 £0.031 0.137 £ 0 020 

I4 n 0.170 £ 0 040 0.182 £ 0 028 

We chose a consistency statistic for helium that was different than for ele- 
ments with charges Z > 6. Ail of the helium measurements were made with the 
three parameter analysis method. We therefore chose the difference of the two 
mass estimators, 

lint/ = Uj-i.f - . (3.3.4 3) 

where j is the detector number, as a consistency estimator. For the entire HIST 

helium data set, histograms of were accumulated for each RANGE, and At 

to determine the M+jf peak location and width. Helium events were accepted 


i 


J 
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into the high resolution data set If M+// for each event w.u less than i»o stan- 
dard deviations from the peak location. The limits adopted are listed in table 
3.3.43. 

Thbto 3.3.43 - 
Hehum Requirements 
HANGS Limtta(anu) 

1 -0 027 t 0.504 

2 -0.006 ± 0 445 

3 0 023 ± 0.440 

4 -0.001 ± 0.454 

5 -0.014 ± 0 364 


3*3.5. Maas Resolute 

Events that passed the selection criterion of section 3.3. 1-4 were summed 
into a histogram for each mass estimator. The histograms were fit usin,; Gaus- 
sian peak shapes, and the peak location and ay were determined for eat h mass 
estimator. For each RANGE a weighted sum of the mass ejtimators was 1 ormed, 
using for weights ( oy )~ *. as determined in the fits of the individual mass estima- 
tors. Table 3.3.5. 1 lists the mass resolution obtained by element and RANGE, 
where Wj is the weighted sum for RANGE j. 

The mass resolution that we measured for the HIST flight <iata is 
significantly larger than the nominal resolution calculated in chapter 'i, table 
2.0.5. 1. Let us compare the measurements cf ay for U Bx , U xi , and to the 
calculations of ay for R/ L = 2 and RANGES 1-3 in table 2.8 5 1. The experimen- 
tal measurement of ay for helium is 35% to 50% larger, ay for carbon i i about 
65% larger, and ay for magnesium is 40% to 75% larger, than was calc ulated 
The reasons for this discrepancy are not known at present. However, the preci- 
sion of the main results of this work, involving the abur dances of isotopes *He. 
IS C. **Ne. “Mg. and “Mg. Is mainly limited by the statistical fluctuation!: In the 
number of particles observed in the 78 266 solar flare. 
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Table 3.16.1 - 

HPT Mae ■ Ree olullon 


Mast 

Helium 

Max 

0.152 ±0.003 

Mb x 

0 169 ±0 002 

*x 

0.127 ±0.002 

Mxu 

0 .136 ±0 003 

Mm 

0 .170 ±0 003 

M m 

— 


0.122 ± 0.002 

Mm 

0 .139 ± 0 003 

M ia 

0 173 ±0 003 

Mm 

— 


0.124 ± 0 003 

M a 4 

0.113 ± 0 003 

M u 

0 108 ±0 006 

M\4 

— 

*4 

0.099 ± 0 003 

Mm 

0 094 ± 0 003 

Mu 

0.155 ±0 005 

** 

0.088 ± 0.003 

Table 3.0.5. 1 - 


(contlnued)- 


Maaa 

Oxygen 

Mai 

0 307 ± 0 006 

Max 

0 293 ± 0 000 

"i 

0.240 ± 0 005 

My, 

0.261 ± 0 006 

M n 

0 300 ± 0 007 

M At 

0 2 P 7 ± 0 008 

*t 

0.178 ± 0 005 

M a 

0 214 ± 0 005 

M, a 

0.270 ± 0 007 

Mat 

0.344 ± 0 010 

*9 

0.170 ± 0 005 

Mu 

C .152 ± 0 012 

Mu 

0.285 ± 0 024 

Mu 

0.360 ± 040 

*4 

0 137 ± 0 012 


Carbon Nttrogen 

0.225 ± 0 006 0 216 ± 0 017 

0 228 ± 0.007 0 277 ± 0 025 

0.194 ±0.006 0.101 ±0.017 

0 203 ± 0.006 0.208 ±0.018 
0 256 ±0 000 0 211 ±0.018 

0.253 ± 0.009 0.272 ± 0 028 

0.159 ±0.006 0 162 ±0 013 

0.191 ±0.007 0 105 ±0 018 

0 208 ± 0 008 0 232 ± 0 022 

0.303 ±0011 0 276 ± 0 033 

0.159 ±0 006 0.137 ±0 013 

0.148 ±0 010 0.130 ± 0 006 

0.227 ± 0 018 0 273 ± 0 012 

0.322 ± 0 024 0 373 ± 0 020 

0.120 ±0 009 0 146 ± 0 008 


Neon Magnesium 

0 377 ± 0.031 
0.302 ± 0.024 
0.268 ± 0.020 

0 204 ± 0 024 0 310 ± 0.032 

0 288 ± 0 032 0 401 ± 0 052 

0.350 ± 0 046 0 461 ±0 131 

0.170 ± 0 019 0 232 ±0.024 

0 213 t 0 020 0 230 ± 0 029 

0 200 ± 0 027 0 282 ± 0 042 

— 0.285 ± 0 057 

0.203 ± 0 020 0 211 ± 0 027 

0.224 ± 0 041 0 233 ± 0 019 

0 399 ± 0 069 0 363 ± 0 057 

0 232 ± 0 048 0.194 ± 0 018 
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3l4. A Critical Evaluation of tho Itch Resolution Data Sat 

The high resolution date eel Is to be used for the purposes of computing 
the ratio of abundances of two Isotopes of the same element In all the cases to 
be considered here, especially for the elements He. C, N. and 0, one Isotope is 
very much leas abundant than the other. Much of the difficulty In estimating 
the Isotope abundance ratios arises from the necessity of estimating the shapes 
of the distribution functions for the mass estimators, due to a single isotope. 

Unfortunately, HIST was not exposed to pure beams of single Isotopes for 
calibration purposes, and the flight data must be used both to calibrate HIST 
and to measure isotope abundances. We will make plausible inferences about 
the shape of the various mass estimator distribution functions in order to 
obtain isotope abundance ratio measurements 

Carbon 

The measurement of the ,8 C/ ,f C abundance ratio is the most difficult 
measurement attempted here Figure 3 4 1 shows a plot of Mb\ vs U41 for 
RANGE 1 carbon. There is a background source for events that have high values 
of Mb\. which may interfere with the clump of events at (13,13) that are prob- 
ably due to the Isotope l8 C. The events with high U B \ values are possibly caused 
by M2 retriggering, a phenomenon that is discussed in Spalding (1983). 
Retriggenng is a type of error that can occur in the analog to digital conversion 
process. If the output of an ADC is equal to ga.in'uolta.gt 4- offset, then 
retriggering error gives a digital output equal to gaxn*uolta.ga + 2^>//sef, 
where o//sef is about « 30 channels. 

Figures 3.4 2a-e show the weighted sum mass histograms for carbon events 
in the time Interval 78:286:10 to 78:272:00 (the 78:286 flare), for each RANGE 
and for the sum of all RANGES We will take a very simple approach to finding 
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the ,a C/ lt C abundance ratio We will count the number of events in the inter- 
vale 12.5 * W t * 13.5 and 11.5* * 12.5 and assign them to " ,a C" and 

respectively. Then we will try to correct for spillover from the '*C peak only, as 
spillover from the IS C peak Into the (t C peak makes a statistically insignificant 
contribution to the IS C/ ,B C abundance ratio. 

In RANGES 1 and 2 it appears that a smooth continuation of the l, C peaks 
will show that events with 12.5 * Wj * 12.0 are most likely ll C events. Thus we 
excluded events with 11.5 * Wj * 11.6 from the IB C events and 12.5 * Wj * 12.6 
from the **C events. For the remainder of the events In the mass Interval 
12.6 * Wj * 13.5, we made plausible assumptions about the minimum and max- 
imum number of 18 C events In the mass Interval, and listed those assumptions in 
table 3.4.1. For example, in RANGE 1, a smooth continuation of the ta C peak 
might also Include up to 2 additional events in the Inter val 12.6 * Wi * 13.5. 

The number of events around mass 11 Is comparable to or even greater 
than the number at mass 13. The two situations are not the same because, as 
seen in the mass vs. mass cross plots, energy loss mechanisms like those that 
operate In double strip events, or high radius events, caut*. the mass estimators 
to be lower, increasing the low mass background. Thus it makes sense to con- 
sider the background separately above and below the main isotope peak. 


Table 3.4.1 - 


RANGE 

1 

2 

3 

4 


Carbon Haas Assignments 


nmx ,a C nXn ia C eat >a C 

13 11 12 

3 2 2.5 

6 4 5 

2 1 1.5 

21 


l, C aampUng interval 

Total C 

u.0* w, * 12.5 

784 

11.8* W 8 * 12.5 

501 

11.5* r s * 12.5 

535 

11.5* W 4 * 12.5 

151 


total Rl-4 


24 


16 


2031 


ORIGINAL PAG£ 13 
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figure 3.4.1 - A plot of ifa rt. U Al for tho time interval 78 206:10 to 
78:272 00, for RANGE 1 carbon events. Tho man comlotoncy requirement* 
aro shown as dashed linos. 
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flfura 3.4.8ft-b - Mui hlitograma of thi wilfhtid a urn man laUmator for 
carbon, for thi tuna intinra) 78:266:10 to 78 272 00 n*uri 3 4 2a H for 
RANGE 1 and Ofuri 3 4.2b ti for RANGE 2 
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figure 3.4.2o-d - Mata histogram of tha walghtad aum mas* aatlmtor for 
carbon, for tha time interval 78:266 10 to 78:272 00. rigura 3.4.2c la for 
RANGE 3 and figure 3 4 2d la for RANGE 4. 
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M|uw 3.4.2* - A man hiitogram of th* w*kght*d Bum mos •■timator for 
carbon, for RANGE* 1*4 and for th* Um* InUrval 78:266:10 to 78 272 00 
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The RANGE 1 cross plot for oxygon In the 76:266 flare la shown In figure 
3.4.3. For ,a O tha background situation la much battar than for >S C, tho >a O 
araa la furthar from tha sourcs of tha high mass U§\ background to urea Flg- 
ura 3.4.4 shows tha weighted sum mass histogram for RANGES 1-4. 

For an upper Unit to tho ,T 0/ 1 *0 abundance ratio, wo noticed that the 
region 16.5 a Wj a 17 0 has 56 events and the region 17.0 a a 17.5 has only 
one event. Then to measure tho ,T 0/ ‘*0 ratio, we used tha Interval 
17.0 a Wj a 17.5 for tha ,T 0 event sample and the interval 16.0 16.5 for 

tha 1# 0 event sample. 

For the '*0/ '*0 ratio, we counted ,# 0 events In the Interval 
15.5 a Wj a 16.5. and ia 0 events In tha Interval 17 5 a Wj a 16.5. For a back- 
ground estimate, there Is one event In tha Interval 10.5 a Wj a 19.5. so we 
Inferred a maximum number of 6. a minimum number of 5 events, and a "best” 
estimate of 5 5 events due to ,a 0 in the Interval 17.5 a Wj a 18.5. 

Table 3A2- 


Qxygen Haas Assignments 


RANGE 

,? 0 

'•o 

lt 0 sampling 

interval 

>•0 

1-4 

1 

... 

16 0 a 1 

a 16.5 

1440 

1-4 

— 

5-6 

15.5 a 1 

l a 16 5 

3142 


-77- 


ORIGINAL 

0F POOR 


I'-'AG.T 13 
Quality 


figure 3.4.3 - A U§ x v* M dl cron plot for RANGE 1 oxygon event* in the time 
interval 78:266:10 to 78:272:00 The dashed lines show the mate con- 
sistency requirements adopted. 
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figura 3.4.4 - A histogram of tha waightad turn mass estimator for oxygen In 
RANGES 1-4 for tha time Interval 78:268 10 to 78:272 00. 
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Mtnftn 

The nitrogen weighted sum histogram for tha 78 266 flare for RANGES 1*4 is 
shown In figure 3.4.3. A smooth continuation of tha U N peak may result In up to 
two events in the Interval 14.5 * Wj * 15.5 due to U N. giving a maximum number 
of 4. a minimum number of 2, and a best estimate of 3 ,B N events There are 360 
events in the interval 13.5 « * 14 5. 
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figure 3.4.5 - A mass histogram of the weighted sum mast estimator for 
nitrogen in the time interval 70:268:10 to 78:272:00. for RANGES 1-4. 
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The neon measurements for the 78:266 flare are shown in Azures 3 4 6a-b. 
To calculate the "Ne/ "Ne abundance ratio we uoed data from RANGES 1*4. The 
"Ne measurements are not influenced as much as carbon and oxygen, by the 
high mass shape of the main Isotope peak, and we estimate that the back- 
ground at "Ne due to "Ne is negligible. To calculate the B, Ne/"Ne abundance 
ratio we will use only RANGES 2-4, to avoid making any guesses about the RANGE 
1 “Ne peak shape. 

Thble 3.4.3 - 


RANGE 

"Ne 

"Ne 

"Ne sampling interval 

"Ne 

1 

0-6 

12 

19.5 £ W t £ 20.5 

103 

2 

0 

12 

19.5 £ IT, £ c0 5 

62 

3 

0 

7 

19.5 £ JT S £ 20.5 

77 

4 

0 

0 

19.5 £ W 4 < 20.5 

6 

total R2-4 

0 

19 

19 5 £ £ 20 5 

145 

total Rl-4 

0-6 

31 

19 5 £ £ 20.5 

246 



0fur* 3.4.8e-b - Neon weighted a urn mail estimators for the time Interval 
78:286:10 to 78:272:00 Figure a la for RANGE 1 and figure b Is the sum of 
RANGES 2*4. 
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A histogram of the weighted aum mass estimator for magnesium la shown In 
figure 3 4.7. To estimate the “Mg/ **Mg and “Mg/ ,4 Mg abundance ratios we will 
assume that the single Isotope peak shape Is a Gaussian. “Mg and “Mg are 
abundant enough that the shape of the tails of the mass estimator distribution 
function Is not very Important. In section 3 5 we will describe the method used 
to obtain rallies and uncertainties fdr the Mg abundance ratios. 


figure 3.4.7 - Magnesium weighted sum mass estimator In the time interval 
78:268:10 to 78:272:00. for RANGES 2-4. 
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HtUum 

For helium we analyzed the data in two parte, RANGE 1 and RANGEe 3-5. 
The data for the weighted sura results for RANGES 2-5, In the time interval 
78:206:10 to 78:272:00, are shown in figure 3.4 8b. We chose mass sampling 
intervals of 2.5 < Wj * 3.4 for *He and 3.5 * Wj * 4.4 for 4 He because we felt that 
It was likely that events in the mass interval 3.4* Wj * 3.5 were due to 4 He. For 
RANGES 2-5 we found no 9 He events and 1020 4 He events. 

The weighted sum results for helium RANGE 1 events in the same time 
interval are shown In figure 3.4.8a. Figure 3.4.9 is a plot of weighted sum RANGE 
1 data for the *He-rich time periods 78:266:09 to 76:297:06 and 78:307 00 to 
78:306:00, for comparison and to establish the accuracy of the mass scale. 
Using the same mass sampling intervals as for RANGES 2-5, we find 4 *He events 
and 512 4 He events in the data from the 78:268 flare. 

The RANGE 1 data set is of poorer quality than the RANGES 2-5 He data set. 
We would like to consider the possibility that the events in the 9 He RANGE 1 
mass sampling interval are caused by 4 He particles. If the 4 He particles may 
have experienced simultaneous "signal defects" in the measurements of the 
energies deposited in detectors Ml and M2, they might have mimicked 9 He par- 
ticles with no "signal defects". Figure 3.4.10 shows a cross plot of U Bx vs. M A , 
for the RANGE 1 data, with the consistency requirements drawn in as dashed 
lines. There is a group of three events near the position (3,3), which may be the 
result of simultaneous signal losses in detectors Ml and M2. 

If we assume that the probability of a signal defect in M2 is independent of 
the probability of a signal defect in Ml, a rough estimate of the background due 
to 4 He around the position (3.3) can be made If we restrict U Bx to the interval 
3-5 * U Bl * 4.5, then the ratio of events with 3.5 £ A/,u < 4.5 to events with 
2.5 * U Al * 3.5 is 0.075 t .012. For the background estimate we restrict U Bx to 
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2.5 * M§ i < 3 5. There la a small correlation between a signal defect in ^elector 
M2 and the mass estimator Mai- We expect the main peak of the Max histogram 
to shift slightly, under the requirement for a signal defect in M2, which is 
represented by the equation 2.5 < M§\ * 3.5. Then there are 40 events with 
2.5 < Mb\ < 3.5 and 3.4 * Mai * 4.4. with a mean mass of 3.89 ± .03. The number 
of events predicted in the region 2.4 « Jtfji < 3.4 and 2.5 « Mg x < 3.5 is then 
3.0 (43.0. —1.7), which is consistent with the observed number of events. Thus 
we will evaluate helium observations as placing an upper limit on the amount of 
®He observed. Table 3.4.4 displays the observed numbers of He events. 


Tfcble 3.4.4 - 


RANGE 

*He 

4 He 

4 He sampling interval 

1 

4 

512 

3.5 < W, < 4.4 

2-5 

0 

1020 

3.5 * W, < 4.4 


ORIGINAL PAGE IS 
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figure 3.4.8*-b - Helium weighted sum man estimator in the time interval 
76:200.10 to 76:272:00. Figure 3.4.0a Is data from RANGE 1. and figure 
3.4 8b Is data from RANGEs 2-5 
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Acur* 3.4.9 - Helium weighted sum oua estimator for RANGE 1 . The tim* 
intervals are from 78 298 09 to 78:297 09 and 307:00 to 308 00, and include 
two *He-rtch flares. 
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Ofura 3.4.10 - Jfa n U Ai for RANGE 1 helium •▼•nta in the time interval 
76266:10 to 76:272 00 Tho daahed Unot show tho rnaaa consistency 
roqulromonta adopted. 
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Tor each element. the result* of tho mui measurements wars expressed as 
the ratio of tho obundonco of o rarer Isotope to that of the most abundant Iso- 
tope for that element. For nitrogen, the result labeled as "'•N/ U N" is the ratio 
of the observed abundance of 1V N to the abundance of **N. The maximum likeli- 
hood technique was used to determine the most likely value for the abundance 

ratios, and to calculate the statistical uncertainties In the abundance ratios. 

* 

For the element magnesium, the weighted sum histogram for each RANGE 
had been fit by the chi square technique, with Gaussian peak shapes. This 
determined the rms. width, o«, and the location of tho mass peaks. Fixing 
and the location of the mass peaks for each RANGE, the likelihood function had 
two free parameters, a * >B Mg/**llg and fi = **Mg/ M Mg The free parameters 
were allowed to vary to find the maximum likelihood solution, fitting all RANGE* 
simultaneously. The statistical uncertainties In the two parameters were found 
by numerically Integrating the likelihood function Let 

f f L{afi) dad ft 

U (*) = ^ • (5 8 . 1 ) 

£ £ L( aft) dadfl 

Then the confidence Interval for the parameter a at a level of confidence C is 
[ Vi.V*]"k«™ 

(3.5.2) 

and 

l-( 1 f £ ) s /.(Vt) • (35 3) 

The probability that a will be found within the confidence Interval Is C, If the 
priori" probability distribution for a Is uniform on the real axis. 


For magnesium the form of the likelihood function chosen was 
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ft *(*.«*) (3 5.4) 

with 

^ “ rS* ] O S O) 

where 


a+A«-l 


a+0-M 


a+0+1 


(35.6) 


and 


a * ®lif / ^Mg ; fi * "Mg/ ,4 Mg ; 

i/y * man o/ iratop* y , 

M * umightrd rum man rrtimator for particlr i . 

For carbon, tha uncertainties In tha abundanca ratios were dominated by 
uncartainUea in tha shape of the man peak for ,8 C. Wo than made plausible 
assumptions about the mass peak shapa, as detailed In section 34, and 
obtained maximum and minimum abundances for tha iso cope ,a C. The likeli- 
hood function chosen was 


Ua.0) = (A lt )" Uw)*“ Mu)*” . (3 8 7) 

a * ,8 C/ **C ; p * ,4 C/ **C ; 

; Al%s 7&T : x ‘ 4= ^4^T : 

where N Is the total number of events observed and n, 3 la the number of l8 C 
events. The upper limit to the abundance ratio ,8 C/ 18 C was obtained by solving 
equation 3.5.3 with the above likelihood function, with n l3 set equal to the max- 
imum number (24) of t8 C events from section 3.4. The lower limit to the abun- 
dance ratio was found by set f, ng n 13 to the minimum number (16) of 18 C events, 
and solving equation 3.5 2. We used the mean of the maximum and minimum 
number of ,9 C events as our best estimate of the I9 C abundance. 

The elements nitrogen and helium used a likelihood function of the type 
3.5.7, but with only one varl&ble. as these elements have only two long-lived iso- 
topes. 
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Oxygen 1B 0/ ,a O and ,7 0/ ia O abundance ratios and uncertainties were cal- 
culated separately, with one variable likelihood functions of the type 3.5.7 and 
using the data in table 3.4.2. 

The neon zl Ne/ zo Ne abundance ratio was estimated using the data from 
RANGES 2-4, with a likelihood function similar to equation 3.5.7. The ^Ne/^Ne 
abundance ratio and uncertainties were calculated with the data from RANGES 
1-4 and a likelihood function similar to equation 3 5.7. 

The results of the confidence interval calculations for all elements 
analyzed appear in table 4.1. 


3.6. Isotope Sampling Corrections 

Different isotopes of the same charge, that stop in a given RANGE, sample 
energy intervals that differ slightly. If we approximate the range energy rela- 
tion by a power law. we then find that the isotopes A! and <4 Z . which each stop in 
a distance R,, and have energies per nucleon and X^z. are related by the 
equation 

Xa\ _ 

%A2 ^1 

Thus the heavier isotope has a lower energy per nucleon to travel the same dis- 
tance. This effect introduces a bias into any isotope abundance ratio measure- 
ment obtained by integrating measurements made over one RANGE. 


(3.01) 


We can compute the bias in the isotope abundance ratio introduced by 
integrating the observations by RANGE, ii we assume that the isotopes in the 
abundance ratio have the same spectral "shape*’. Suppose that the flux of iso- 
tope j, per cm 2 — sec — sr - Afe V/mic , at an energy/nuc, e, is 

0. (3 8 1) 


and that the flux for isotope i is 
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=*/(•)• (3 8 2) 

For element X. the ratio of the abundance of isotope i to the abundance of iso- 
tope j, is 


*L = *L 

Xi Ki 


(30.3) 


and the spectra of isotopes i and j are considered to have the same f pectral 
shape as a function of energy/nuc. 


Let < Xj> be the number of events of isotope j. measured in a given RANGE, 
an interval in pathlength from R(L) to R(U). Let e^(U) be the energy/nuc, such 
that isotope j is stopped by a thickness R(U) of silicon- Th?n 

<Xj>= f Kj /(■)*« zKjC, (3 6 4) 

•,(t) 


and 


Kj__ <Xj> q 
Xi ~ < Xi> q 


(3.0 5) 


Thus we must multiply the observed ratio of < Xj >/ < X i > in a RANGE by the 
correction factor q/ q. Suppose that OT Ne and ^Ne heve the same : pectral 

shape as a function of energy/nuc, and that /(e) = (~) 3 If R(L) is equt l to the 

range of ®”Ne at 25 MeV/nuc, and R(U) = <*>, then 


X a < Z2 Ne> 


0 897 


( 3.0 6 ) 


Xgg < w Ne> 

We can measure the spectra of the more abundant iso’ ope as a function of 
energy/nuc, AJ/(e). If we assume that the 3r ectra of the less abund int iso- 
tope has the same shape as a function of energy/nuc, we can compute the 
correction factor Cj/ C\. That assumption is equivalent tc the assumption that 
the isotopes’ spectra have the same shape as some arbitrary function of the 
isotopes' velocity 


The diffusive propagation of solar energetic particles from Sun to Ilarth in 
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the interplanetary medium is thought to depend on the particles' rigidity, or 
momentum per charge Thus two isotopes of the same element night have 
spectra which were the same shape as a function of rigidity We can also com- 
pute the corrections to the isotope abundance ratios for this case If the two 
Isotopes carry the same number of electrons, then they also have the same 
spectral shape as a function of their total momentum. 

Suppose that the (lux of isotope L per cm* — sec - sr - MeV/ c is 

^- = *>y(P) (36.7) 

and the flux of isotope i is 


fjr = *i9(P). (3.6.8) 

Let Pj (U) be the total momentum, such that isotope j is stopped by a thickness 
R(U) of silicon. Then 


p, cm 

< Xj> = f A)c/(P)dP . 
P, CL) 


If p is equal to the momentum per nucleon, pj = 


V 


For isotope i. 


P){ u) 

< Xj> — Aj J Kj g{Ajp) dp s Kj A, Dj , 

J»yO) 


(3 6 9) 


(3 0 10) 


i\ on 

<X t >= f A;p(P)dP (3 0 11) 

1\ (L) 

*(m 

= A f K i g{A i p)dp 
*04 

(V-V*(t0 

- A i f A \g{A,p)dp 

IV 4) F|W 

= 4 ■ 

Thus we have 


X} _ < j 

* <*> iDyy 


(3 0 12) 
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We can measure the spectra of isotope j, as a function of momentum per 
nucleon. 


di 


dp 


= 44 3 (4p) 


(3.8.13) 


dj, 

If ire assume that — has the same spectral shape as a function of total 
a P 

a i 

momentum or rigidity, then we can calculate and correct for the RANGE 


sampling effect, as follows. Let 


rfj, 

X * - *• * •*> • 


and 


(At/ Aj) p^J) 

r 

(V^)PiCU dp 


f, - ^ * ■ *! *> 


then 


(3.6.14) 


(3.8.15) 


Wj Dt 




^ J 


(3 8 18) 

Suppose that ^Ne and ^e have the same spectral shape as l function of total 
momentum, and that 


dJ gQ _ f l_\s de(p) 
dp V e ‘ dp 

For R{L) equal to the range of at 25 MeV/nuc. and /?( U) = <■>. 


(3 6 17) 


1T~ = tSt^ 1308 - ( 3810 ) 

* 20 < Ne> 

Note that these two assumptions about the spectral shape have produced 
correction factors that operate in different directions upon the data 


The isotope abundance correction factors were calculated for each range 
and isotope ratio, using both of the above assumptions about the functional 
dependence of the spectra The spectral shape chosen for elements C N. 0. Ne, 
and Mg. for the most abundant Lsotope. was the six element spectral fit of sec- 
tion 3 9. For helium, we used the helium spectral fit, also from section 3.8 The 
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integrals of equation 3.8 2 were done numerically. 

To obtain the isotope abundance factors for measurements that integrated 
data over several RANGES, we computed a weighted mean of the correction fac- 
tors for each RANGE, using for: weights the number of events of the given ele- 
ment observed in each RANGE 


The isotope abundance correction factors are listed in table 3.6.1. 
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Tfeble 3.6. 1 - 


Ratio 

RANGE 

R/nuc 

functional dependence 

Momentum 
functional dependence 

,S C/ ia C 

1 

0 942 

1.131 


2 

0 910 

1.243 


3 

0.883 

1.371 


4 

0.824 

1 673 


1-4 

0.911 

1 288 

,4 C/ ,8 C 

1 

• 0 892 

1 292 


2 

0.852 

1.560 


3 

0.790 

1.871 


4 

0 894 

2848 


1-4 

0 839 

1.S30 

i5 N/ u n 

1 

0.935 

1.120 


2 

0.919 

1.232 


3 

0 891 

1.327 


4 

0.S34 

1.625 


1-4 

0.913 

1.230 

,7 0/ lfl 0 

1 

0 950 

1.121 


2 

0.919 

1.214 


3 

0 899 

1.299 


4 

0.836 

1.537 


1-4 

0 928 

1.195 

ibq/ ie Q 

1 

0.907 

1.264 


2 

0.852 

1.470 


3 

0.813 

1.695 


4 

0.713 

2.178 


1-4 

0 087 

1.438 

z, Ne/ eo Ne 

1 

0 948 

1.116 


2 

0 929 

1.180 


3 

0 908 

1.263 


4 

0 079 

1.319 


2-4 

0 918 

1.228 

“Ne/ eo Ne 

1 

0.902 

1.252 


2 

0 886 

1.391 


3 

0 914 

1.605 


4 

0.773 

1 629 


1-4 

0 089 

1.400 
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Tfeble 3.6.1 - 
(coni.) 



K/nue 

Momentum 

Ratio RANGE 

functional dependence 

functional dependence 

Mg/^Mg 2 

0.935 

1.158 

3 

0.914 

1.241 

4 

0.921 

1.182 

1-4 

0.923 

1.202 

Mg/ 24 Mg 2 

0.877 

1.342 

3 

0.840 

1.540 

4 

• 0.848 

1.335 

1-4 

0.858 

1.444 

3 He/ 4 He 1 

1.247 

0.873 

2 

1.310 

0.597 

3 

1.404 

0.504 

4 

1.470 

0.360 

5 

1.874 

0.259 

1-5 

1.348 

0.568 

2-5 

1.398 

0.510 


3.7. Low Resolution Data Set 

We defined a low resolution data set for the purpose of measuring the 
energy spectra of various elements measured by HIST. We relaxed the require- 
ments of the high resolution data set to obtain more events and more coverage 
of lower energies (primarily through the acceptance of RANGE 0 events) 

Element identification was through two parameter analysis In RANGES 0 
and 1 we required that the mass computed using the '.wo deepest detectors 
triggered be within 3 amu of the principal isotope of the element. In RANGE 1, 
carbon events thus had to have 9 -s U Bl « 15 to be included in the low resolu- 
tion data set. For RANGES 2-5, the two parameter mass had to be equal to the 
principal mass plus 3 or minus 2 amu Single strip triggers were required in 
matrix detectors used for the mas3 estimators (Ml and M2 for RANGE 0. and M2 
for RANGE 1). No energy limits were placed on El, the energy deposited in the 
stopping detector The radius limits were the nominal detector radii from the 
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detector specifications 

For the element helium, the mass limits were 1 and 6 amu for all RANGES 
For helium in RANGE 1. retriggering in detector M2 (see Spalding 1 9t 3) made 
some protons have mass measurement errors which made the protons indistin- 
guishable from RANGE 1 helium ions, using the Mgi mass estimator The protons 
and helium ions were then separated by requiring that 10s Ua\ ^ 6 0 .Uso hold 
for RANGE 1 helium 

• " 

For iron we could not require single strip triggers in the matrix detectors 
used for mass estimator:? and t ' able to use the hodoscope efficiency calibra- 
tion of appendix A. The hodoscope efficiency calibration was not extended up to 
the energies characteristic of iron ic as stopping in a matrix detector For- 
tunately, signal defects in the matrix detectors for iron events do not iDwer the 
mass estimators enough to make it impossible to differentiate between iron and 
calcium. An iron event in RANGE 0 or RANGE 1 with a signal defect in VA or M2, 
still lies above the "track*' for calcium events We then chose mass corsistency 
requirements Tor "iron" events that excluded calcium events but sumrcied over 
the charge interval 21 ^ Z ^ 28 If the SEP abundances were smilar to 
Cameron’s 1981 abundance table in this charge region, then this procedure 
would result in an 8% overestimate of the SEP iron abundance, much less than 
our experimental errors. 

3>8. Element Spectra 

Using data from the low resolution data set we computed the spectra of 
elements He,C,N,0,Ne,Mg,Si, and Fe During the highest rate time periods, a 
number of anomalies in the operation of HIST occurred, which are desmbed in 
Spalding 1983 HIST was operating at counting rates very near to the maximum 
rates at which it was designed to function correctly, and the causes of ill of the 
high rate anomalies are not completely understood However, we bebeve that 
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we have been able to interpret the data correctly and extract the physically 
interesting information from the data 

We compute fluxes for time periods that are an integral number of 3 hour 
time periods Let the 3 hour time periods be labeled by 1. We compute a flux, 
F^, of a given element k, in an energy Interval i. and for a time interval j, con- 
sisting of 3 hour time intervals from l - lmm(;') to f = fmox(j'), 

T7 Hu* 

r* - r . • C3.B-D 

/ f Dsnom^ 

l = 

with 


Lkmom^ = [^{An^CORRj ^(RLTi^EROj L ) (3.0.2) 

where 

Nut = number of events, of element k. in the low resolu.ion data set, 
in the energy interval i, and the 3 hour time interval l 

A Ei = the width of energy interval i, in MeV/nuc. 

Af)/ = area solid angle geometry factor for RANGE I. 

(calculated in Mewaldt (I960)) 

/(i) = RANGE of energy interval i. 

RLTi = rate live time is seconds, 3 hour time interval l. 

For H1Z events, the event read out efficiency is given by, 

_ number of HJZ events of RANGE I 
1 number of HJZ rate counts of RANGE I 

LQZ rate events are not accumulated continuously by RANGE, so for LOZ events, 

/rpn - number of LOZ events 
^ number of LOZ rule counts 

Suppose that the events in interval ilk are labeled by q, then the hodo- 
scope sampling correction for RANGE 0, CORK/ L k , is 

CORR llk = -±- £ [(1 -/Wtfl.£in(g)))(l (»)))]"' 


-°-r P 00«O^TV 

For RANGE 1. 

CORR,t> = “ ? (i - 

and for RANGES 2-5, 

CORR Iik - 1.0 

For "iron” events CORR = 1.0. 

Here, ,E Mx {q)) is the probability of a double stnp trigger in detec- 

i * 

tor Ml. for a particle q that deprshed an energy Eh 1 ( 7 ) in detector Ml. The 
P ‘toil function is described in appe/ tix A. 

3.9. Fits to the Element Spectra 

To determine a spectral fora. /(*). with which to compute isotope abun- 
dance ratio corrections, we averaged data over the 78 206 flare For the time 
period 78:266:10 to 78:272 00 and for elements C, N, 0, Ne, Mg, and Si. average 
fluxes were computed in two energy bins per RANGE, for RANGES 0 to 4 All six 
sets of data were fit to a common spectral form by the chisquare technique 

The average flux for element j, between v ^ and energy /nuc is F^j 

The spectral form, /(z). was defined by ten free parameters The parameters 
were the value of / , |y t {. at 10 points, }z,j. equally spaced in log energy from 4 
to 80 MeV/nuc / was interpolated log-linearly between the energy points, 
if 

k = integer part 0 } j log (z )- (09 (4) J/ log del . (3 9 1) 

with 

logdel = ( log (60) - log (4))/ 9 . (3 9 2) 

then 

!»,(/(,))= Egill -Vwy - ‘o jilll ;tsg( V..,) - log ( v«)l 1333) 


+ y*) 
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The abundance of each element was normalized to silicon The flt flux of 
element j, at energy /nuc 2 . was given by P ,/(*). with ^ - 10 The 15 free 
parameters. (y»j and \p } \, were found by minimizing the function. 

t 2 \fii9* . (3 9 4) 

J-C4-1 ' 

where 

*ui j 

9m ~ f /(*) d« / . (3 9 5) 

*» 

and a w is the statistical uncertainty ir the flux F^ 

The data and the resulting flts for ail the elements flt are shown in figures 
3 9.1a-f is plotted as a solid line, P } g^ is plotted as a dotted line, and 

the fluxes A v are plotted as error bars The data support the hypothesis of a 
comriK>n spectral form The main systematic deviations from the flt are the 
higher energy points of the element 3ilicon, wtich show higher fluxes than the 
common spectral form would allow Figure 3 9 2 shows the data points super- 
posed, wo plotted F^/ vs energy / nuc to show the common spectral shape 

For the element helium, the fluxes averaged over the 78 266 flare were flt 
to a spectral form f [p) = A axp[-p/ p 9 ), where p is in units of momentum per 
nucleon, using two free parameters, A and p B We found p 9 = 2b 54 t 0 43 The 
flt and data are shown in flgure 3 9 3 

These spectral flts were used in section 3 6 to compute isotope abundance 
correction factors Fits to the element abundances on a daily basis are 


reported in section 4 3 
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flgurr 3.9.1a-f - A plot of the average spectral flux, Fy, pei 
cm* - sr - sec - MaV/nx ic , for the time interval 70 266 10 to 78 272:00 
The solid lines are the spectral form, pjf(E/n). The dotted lines are the 
fit integrated over the experimental energy bins, Pjihj The experimental 
points. Fy. are plotted with error bars 
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figure 3.9.2 - A plot of the average spectral flux for elements C. N, 0. Ne. Mg. 
and Si. F^/ pj, for the time interval 78:268:10 to 78 272:00. The fluxes of 
elements C, N. 0. Ne. and Mg are normalized to silicon. 
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figure 3.9.3 - A plot of the average spectral flux per 

cm* - sr - sec - (fJeV/ c)/ nuc . vs momentum per nucleon for helium, in 
the time interval 78 260 10 to 78:272:00 The plot symbols are the same as 
for figure 3.9 la-f 
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Chapter 4 - Observations 


4.1. Owview 

The HIST instrument was operational from 1978 22b UT to 1978 335 UT On 
1978:335 UT a failure occurred in the digital electronics (see Mewaldt 1980) 
that read out events into the telemetry stream The failure made HIST unus- 
able as an isotope spectrometer after 1978 335, using the methods of this 
thesis. Breneman (1980) has analyzed a methodology that allows one to make 
measurements with the degraded HIST instrument, at reduced elemental and 
isotopic resolution, and over more limited energy ranges 

h'igures 4 1.1 a-c show the HIZ rate, corrected for live lime, in RANOKs 0. 2, 
and 4, for the time interval 78 225 to 78 vl35 The HIZ rales primarily count C 
and 0 nuclei and table 4.1.1 lists the energy intervals for the HIZ rates, assum- 
ing that all of the particles incident are oxygen The carbon energy intervals 

would be about 15* lower 

Table 4.1.1 - 


(MeV/nuc) 

KANCK Oxygen ene rgjr interval 

0 48- 69 

1 6 9 - 10 ! 

2 10 1 - 14 2 

3 14 2- 24 3 

4 24 3 - 45 3 

In figures 4 1 la-c there are several sharp increases in the counting rate, 
or particle everts, superposed on a relatively constant background There are 
two very small particle events that occur on days 78 278 and 78 283 Two larger 
events occur starting on days 78 260 and 70 314 All of the measurements 
reported in this thesis will be from the particle event that begins on day 78 266 
The event on day 78 314 has too few particles in its high resolution data set to 
allow us to make meaningful conclusions about its isotopic composition Kor 
example, for days 70 314 through 78317, the high resolution data set contains 
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about 30 RANGE 1. 7 RANGE 2, and 5 RANGE 3 carbon events, compared to 
almost 900 RANGE 1. 600 RANGE 2, 600 RANGE 3 carbon events for the days 
70 266 to 76 272 

Figures 4 1 2a-c show the H1Z rate in RANGF^ 0. 2, and 4 for days 70 265 to 
78 275. The event is characterized by a fast rise, an almost flat peak, and an 
exponential decrease The shape of the particle event is approximately similar 
in each range, though the rate of decrease after day 78 268 is faster for the 
higher ranges All of the isotope and element observations reported in chapter 
4 are for the time interval from 78 266 10 to 78 272:00. For that time interval, 
the RANGE 0 H1Z rate is higher than the levels before and after the particle 
event by at least a factor of two 

Slightly before the particle event, a large solar flare of importance 3B in H a 
was observed on the sun. from 78 266 09 44 to 78 266 12:15 (The data in this 
paragraph are from the Solar Geophysical Data Bulletin 1979. March. April ) The 
flaring activity was observed in the McMath plage group 15543, from which 
other smaller H m flares were recorded in the time interval from 78 266 05 to 
78 269 12 From 78 267 02 07 to 267 02 47 a flare of importance IN was seen 
and from 78 267 06 23 to 70 267 06 29 a flare of importance 2B was seen in plage 
group 15543 The Deep River neutron monitor, with a 1 03 GV geomagnetic 
cutoff, recorded the onset of a ground level event at about 78 266 10 30 X-rays 
from 0 5 to 4 A show a sharp increase of a factor of about 600 at 78 266 10 and 
a smooth decline until 78 266 20 Electrons from 1 to 5 MeV increase from 
about 5 x 10 3 /cm 2 sr sec \h*V before the event to a peak of about 
2x 10 M / cm 2 sr sec WeF They remain over 1 x 10 H / cm 2 sr see UeV until 
day 70 208, then decrease approximately exponentially to 
1 x 10 1 / cm* st sec M»F on day 78272 McMath 15543 was located at 35® /V, 
50* W on the sun on day 78 266, and at a solar wind speed at earth of 
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w 350 km/ sec . it was well connected by a Parker spiral magnetic field to the 
Earth. 

In figure 4.1.2, there is a sharp decrease in the particle couiting rates in 
the hour before day 78:208, which may be associated with the arrival of a sud- 
den storm commencement on 78:268:07:18 (Solar Geophysical Data Bulletin 
1979, March, April ). Sudden storm commencements are associated with the 
arrival at Earth of interplanetary shock waves in the solar wind plasma, travel- 
ling from the Sun (Svestka 1976) If the sudden commencement was associated 
with the 3B flare of 78:266 09:44, then the flare shock travelled at an average 
speed of 910 km/ sec to reach the Earth, fairly typical for a flare shock 

All of the above data are consistent with the occurrence of a large solar 


flare on the Sun 
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figure 4.1.1e-c - The HIZ rawe for RANGES 0,2. and 4. corrected for live time, 
is plotted as a function of time. Figures a,b, and c, correspond to RANGES 
0, 2, and 4, respectively. 





- Ill - 


ORIGINAL pace ,s 
0F POOR quality 


figure 4.1.2a-c - The H1Z rate for RANGEs 0,2, and 4, corrected for live time, 
is plotted as a function of time. Figures a,b, and c, correspond to RANGEs 
0, 2, and 4, respectively. 
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4.2. Isotope Observations 

The mass measurements of section 3.4 and the confidence interval calcula- 
tions of section 3 5 were corrected for energy per nucleon sampling effects with 
the factors of table 3 0 1. The results, for the time interval 266 10 to 272 00. are 
shown m table 4.2.1 and plotted in figure 4.2.1. The energy intervals Listed in 
table 4 2 1 were calculated for the most abundant isotope in the respective iso- 
tope abundance ratio We also List isotope abundance ratios computed from a 
compilation of solar system isotope abundances (Cameron 1981) For al! of the 
isotope ratios for which we report a finite abundance ratio, we see no statisti- 
cally significant disagreement with the Cameron compilation. 


I 


Table 4.2. 1 - 
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HOT Isotope Ratio Results 
(corrected for energy /nuc sampling effects) 


Ratio 

Observation 

(MeV/nuc) 
Energy Interval 

(1080) 
Cameron value 

RANGES 

*He/ 4 He ' 

* 0 0028 

4.B -32 5 

0 00018 

RANGES 2-5 

*He/ 4 He 

<0 0175 

3.3 -4 8 

0 00018 

RANGES 1 

,a C/ **C 

0 0095 + 0 0042 
- 0 0029 

0 2 - 38.7 

0.0111 

RANGES 1-4 

,4 C/ **C 

< 0 0014 

8 2 - 38.7 

— 

RANGES 1-4 

1» N/ U N 

0 008 + 0 010 
- 0 005 

7.3 - 42.1 

0.0037 

RANGES 1-4 

,7 0/ "0 

< 0 0021 

6 9 - 45.3 

0 00038 

RANGES 1-4 

18 0 / 1#0 

0 0015 + 0 0011 
-0 0007 

8 9 - 45 3 

0 0020 

RANGES 1-4 

* , Ne/ 20 Ne 

<0014 

113 - 51.2 

0 0030 

RANGES 2-4 

“Ne/ Z0 Ne 

0 109 + 0 028 
- 0 019 

8 2 - 51 2 

0 122 

RANGES 1-4 

23 Mg/ 84 Mg 

0.148 + 0 040 
- 0 028 

12 1 - 55.3 

0.128 

RANGES 2-4 

“Mg/ 24 Mg 

0 148 + 0 043 
- 0 025 

12 1 -55.3 

0.141 

RANGES 2-4 


Isotope Abundance Ratio 
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O^ure 4.2.1 - Isotope abundance results for the 268 flare Plotted as dashed 
lines are isotope ratios calculated from the Cameron (1981) compilation 




- 116 


To derive the spectral corrections to the isotope ratios, we had to make the 
assumption that different isotopes of the same element had the same spectral 
shape. We can now check this assumption. For the isotope ratios ,3 C/ ’*0, 
,B 0/ ,8 0. “Ne/ ,ao Ne. “Mg/^Mg, and the isotope ratio for the tune 

interval 286:10 to 272:00 was computed in each RANGE separately. We plotted 
the abundance ratio results vs. energy per nucleon in figures 4.2 2a-e. If each 
isotope, of a given element, had the same spectral shape, then the abundance 
ratio should stay constant as a function of energy. We see no evidence of a sys- 
tematic variation with energy of the isotope abundance ratios reported here. 
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fl(ura 4.2.2a - Isotope abundance ratios are plotted vs. energy per nucleon. 
Each plotted point corresponds to the results from one RANGE. The hor- 
izontal dashed lines and the arrow correspond, respectively, to the 
confidence limits and most likely values from table 4.2.1. 




I8 0/ I6 0 
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figure 4.2.2b - Isotope abundance ratios ere plotted vs. energy per r ucleon 
Ea^h plotted point corresponds to the results from one RANGE. The hor- 
izontal dashed lines and the arrow correspond, respectlve’y. to the 
confidence limits and mo*,t likely values from table 4.2.1. The dashed 
extensions of the lower error bars for figure 4.2.2b are the "conservative 
error tars discussed in section 1.2. 
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figure 4.2.2c - Isotope abundance ratios are p’otted vs energy per nucleon. 
Each plotted point corresponds to the results from one RANGE. The hor- 
izontal dashed lines and the arrow correspond, respectively, to the 
confidence limits and most likely values from table 4.2.1 
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Another assumption of the above analysis was that the isotope abundance 
ratio is constant as a function of time . ,3 C/ ,l C. 1B 0/ ,a 0, ^Ne/^Ne, ^Mg/^Mg, 
and “Mg/^Mg abundance ratio observations were binned into one day time 
intervals and plotted vs. time, in figures 4.2.3a-e. Isotope abundance correction 
factors, calculated using the methods of section 3.8 for an energy/nuc func- 
tional dependence, were computed on a daily basis from the daily C-Si fits of 
section 4.3. If the abundance of on,e isotope of a given element had a different 
behavior as a function of time than the abundance of another isotope of the 
same element, we would see a systematic variation in the plotted isotope ratio 
as a function of time. If the propagation of the nuclei from sun to earth was 
rigidity dependent as discussed by Scholer at a l. (1978) and Witte at al. (1979). 
then we might see all of the heavier isotopes simultaneously enhanced or 
depleted together at the beginning of the particle event. We see no evidence 
for time variations of the isotope abundance ratios. 
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figure 42.3b - Isotope ratios are plotted as a function of time. The hor- 
izontal dashed lines and the arrow correspond, respectively, to the 
confidence limits and most likely values from table 4.2.1. The dashed 
extensions of the lower error bars for figure 4.2.3b are the "conservative 
error bars discussed in section 4.2. 
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figure 4.2.3c - Isotope ratios are plotted as a function of time. The horizon- 
tal dashed lines and the arrow correspond, respectively, to the confidence 
limi ts and most likely values from table 4 2.1. 
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In figures 4.2.2 and 4.2.3 some of the lower error bars are extended by 
dashed lines. This is to show the effect of the assumption made in section 3 5 
that the abundance ratio is "a priori" equally distributed along the real axis. 
For a very small n umb er of events ( < 3) in the numerator of the abundance 
ratio, this assumption causes the likelihood confidence integral for the abun- 
dance ratio, equation 3.5.2, to be concentrated above the mean value. For 
example, on day 286, 2 l8 0 and 081 ,a 0 nuclei were observed. The value for 
/«(*). where the parameter a is equal to the abundance ratio ,8 0/ ,8 0 and for x 
equal to 2/ 881 is I m {2/ 881) = 0.324. 

A more conservative assumption states that the lower error bar for a 
parameter a, calculated for a given confidence level C, is the value of a such 
that the probability of obtaining more events than are observed is (l-C)/2 
The "conservative” upper error bars calculated by this technique are very 
nearly equal to the upper error bars calculated by the likelihood technique of 
section 3 5 The difference between the two methods for the lower error bars is 
largest for a very small number of events in the numerator of the abundance 
ratio, and we have only plotted the "conservative" error bars for the oxygen 
plots, figures 4.2.2b and 4 2 3b. 

4.3. Hement Abundance Ratios 

The low resolution data set was binned into time intervals of one day and 
average fluxes were computed for each day, for the elements C. N. 0, Ne. Mg. 
and Si. The elements were simultaneously fit to a common spectral form, 
j {p) = exp(-p/ p 0 ). where p is momentum per nucleon and p, is a constant, 
which is the same for all elements Table 4.3.1 lists the values of p t obtained. 
The fits were computed by the chi-square technique of section 3 0. Table 4 3 2 
shows chi-squared per degree of freedom as computed each day The results 
for day 280 and afterwards show that the hypothesis of a common spectral form 
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for elements C, N. 0. Ne. Mg. and Si produces a good flt to the actual spectral 
form Figure 4 3.1 shows the data and the flt for day 268 The error bars for the 
element fluxes were computed using only statistical considerations. Figure 
4.3.2 shows the abundances relative to oxygen obtained by the daily flts as a 
function of time. The abundances, pj, on and after day 268 were weighted using 
the statistical errors obtained in the flt program and averaged to give a best 
estimate of the flare element abundance ratios These values are plotted in 
figure 4.3.2 as dashed lines and listed in table 4.3.3. 
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Table <3.1 - 


Day p, (MeY/c)/nuc 

78:268 29.22 ± 0.30 

78:287 20.25 ±0.19 

78:288 22.44 ± 0.42 

78:289 20.09 ± 0.12 

78:270 19.38 ± 0.96 

78:271 22.13 ± 2.92 


Table 4.3.2 • 


no. tree 


Jmj 

X 1 

no. points 

parameters 

X? 

elements 

286 

2222 

57 

7 

4.44 

C,N,0,Ne,Mg,Si 

267 

341.2 

58 

7 

8.98 

M 

288 

43.0 

52 

7 

0.90 

M 

269 

494 

48 

7 

1.27 

II 

270 

21.0 

35 

8 

0.73 

C,N,0,Ne,Mg 

271 

4.8 

14 

3 

0.44 

C,0 


Table 43.3 - 


Abundance Ratio 

Average 

Comment 

C/O 

0.470 ± 0 018 

average of C-Si fits, 70:288-272 

N/O 

0 125 t 0 009 

" 78:268-271 

Ne/O 

0.132 ± 0 012 

" 78:208-271 

Mg/o 

0 .198 t 0 017 

" 78:288-271 

Si/O 

0 .129 t 0 016 

" 78:268-270 



average of He/0 for 

He/O 

52.8 ±2.1 

RANGES 2-5, 78:269-272 



78:286-272 

Fe/O 

0 137 t 0 030 

RANGE 0 iron, 0 0-10 4 MeV/nuc 

Fe/O 

0 073 ± 0 033 

RANGE 1 iron, 11.8-10 5 l>V/nuc 

Fe/O 

0 034 ± 0.033 

RANGE 2 iron, 16 5-24.3 MeV/nuc 
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figure 4.3.1 - C, N. 0. Ne. Mg, and SI spectra observed on day 260. The fit to 
the data, fj f{p), Is plotted as a straight line. The fit integrated over the 
energy limits of the data points, fij g^. is plotted as a dotted line. The data 
points are plotted as error bars. The flux is in units of 
(cm’jrw {UaV/ c )/ tvuc ) -1 
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I 

figure - Abundances relative to oxygen from the C-Si fit. fij, plotted as 
a function of time. The dashed lines are averages oi the plotted points on 
and after day 200. 

t 
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Helium, especially for days 288 and 287. did not show the same spectral 
shape as the elements carbon through silicon. To obtain a value for the He/O 
abundance ratio, we integrated the oxygen spectra obtained in the C-Si analysis 
over the energy intervals corresponding to each RANGE for helium. Thus we 
obtained an "oxygen" flt flux in each RANGE and for each day that could be 
directly compared to the helium flux. For greater statistical accuracy, pairs of 
RANGES were combined to obtain energy intervals of 2.3 - 4.8 UeV/nuc for 
RANGES 0-1. 4.8 - 11.3 MeV/nuc Tor RANGES 2-3, and 11.3 - 32.5 MeV/nuc for 
RANGES 4-5. The data in the C-Si fits do not extend below about 4.5 MeY/rmc 
for oxygen, so the RANGEs 0-1 He/0 abundance ratios are ba^ed on an extrapo- 
lation of the C-Si flt to lower energies. 

The results obtained are plotted vs. time in figure 4.3.3. The He/0 abun- 
dance ratio varies at the beginning of the particle event and is relatively con- 
stant after day 269 For a best estimate of the He/0 abundance ratio we aver- 
aged the abundance ratios in RANGES 2-5 and for days 269-271 to obtain a value 
of Ha/ O = 52.8 ± 2.1. This value is plotted in figure 4.3 3 as a dashed line and 


listed in table 4.3.3. 
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figure 4.3-3 - He/O abundance ratios are plotted vs time in three energy 
intervals. 
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Iron does not show the same spectral shcpe as the C-Sl elements. Cook 
(1901) has also observed large flares in which the iron spectral shape was 
significantly different than for the C-Si elements. For each RANGE and each 
day. iron fluxes were divided by the "0" fluxes from the C-Si fits, to form Fe/0 
abundance ratios The Fe/0 abundance ratios are plotted vs. tame m figure 
4.3.4. The points for day 288 have a higher value for Fe/0 than does the rest of 
the data set. We averaged the Fe/0 abundance ratio results on and after day 
288 in RANGES 0-2. We And that for RANGE 0. n/0 = 0 137 ± 0 030. RANGE 1 
has Fk/ 0 ~ 0.073 ± 0 033. and RANGE 2 has /k/ O = 0 034 ± 0 033. These 
results are not statistically consistent with a common v^'ue for the Fe/0 abun 
dance ratio The Fe/0 values are listed in table 4 3 3. 
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2-upb 4.3.4 - Fe/O abundance ratios in each RANGE are plotted vs time 
The data points shown as crosses are from RANGE 0. 8 - 10.4 toVVnuc, 
and the points shown as solid circles are from RANGE 1, 
11.8 - 18 5 HaV/nuc. The points shown as Xs are from RANGE 2. 
16.5 - 24.3 Ms V/nuc, and the points shown as open squares are from 
RANGE 3, 24.3 - 43.5 to V/rtuc 



Ds fK OF 



- 138- 


Chapter 5 - Discussion 

5.1. Overview 

In the following section We will review the observations of element abun- 
dances in SEP events. We will show that the 78:266 SEP event has an elemental 
composition that is similar to other large SEP events. We will review the element 
abundance variations and average values observed in large SEP events and 
speculate about plausible consequences for isotope abundance ratios. 

In the third through fifth sections we will review the evidence about the 
Sun’s composition that comes from the observation of solar electromagnetic 
and particle emissions, for the elements He, C, N, 0. Ne, and Mg In the sixth 
section we will review noble gas meteoritic composition and summarize the data 
from sections 3 through 5. in order to compare SEP and solar isotope abun- 
dances. In section seven we will present our conclusions about the spectral 
corrections to the isotope abundance ratios. 

Sx2. Solar Energetic Particle Elemental Composition 

The elemental composition of large SEP events has been observed to vary 
from flare to flare. Since our measurements in this work are limited to one 
large solar flare, we must use other reported solar flare energetic particle 
measurements to establish a context for our observations of the isotopic com- 
position of solar energetic particles. The elemented composition and size of the 
78:266 flare and its location on the sun. well connected to the earth by Parker 
mar.ietic field lines, place it into a well studied class of the largest solar ener- 
getic particle events. The study of a sample of the largest SEP events may pro- 
vide the best estimate of SEP abundances. V.aron et a l. (1980) measured a 
number of element abundance ratios in c largs data set consisting of daily 
averages for the period Oct. 30. 1973 to Dec 1. 1977, in the energy ranges 0 8- 
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1.8 MeV/nuc for H. He, C, and 0. and 1-4.8 MeV/nuc for 0, Ne. Mg. Si. S-Ca. and 
Fe. They find that the abundance ratio variation increases with decreasing flux 
levels. 

Cook (1981) and Cook et a l. (1980) have studied seven large solar ener- 
getic particle events. For four events for which the composition and spectral 
shape are relatively constant with time, Cook (1981) constructed a four flare 
SEP (solar energetic particle) abundance average. SEP(Cook), quoted in table 
5.2.1. Figure 5.2.1. from Cook (1981), shows the element abundances measured 
in each of the four flares used in the four flare average, relative to the 
SEP(Cook) average. The deviation of the element abundances for each indivi- 
dual flare about the four flare SEP(Cook) average, for elements with 6 ^ Z ^ 28, 
was systematic and monotonic. For element abundances normalized to silicon, 
the abundance of the element group C, N, and 0, was observei' to be anticorre- 
lated with the abundance of the element group Fe and Ni. McC uire et a l (1979) 
observed a similar pattern of elemental abundances when comparing SEP event 
abundances with a set of average abundances constructed from SEP events 
with low (iron group)/0 abundances. Meyer (I98!a,b) has attempted to use all 
available data on SEP composition between 1 and 20 MeV/nuc to obtain a "mass 
unbiased baseline". Meyer parametrized the data by the Fe/0 abundance 
observed, and then selected a particular value of Fe/0 so that the abundance 
ratio, SEP/photosphere, for Fe, was approximately equal to the 
SEP/ photosphere abundance ratio for Mg and Si Fe, Mg, and Si all have similar 
first ionization potentials. The average element abundances of 5 event classes, 
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figure 5J2.1 - ( from Cook (1981)) The elemental abundances for four solar 
flares divided by the average abundance for the four flares, SEP(Cook), is 
plotted vs. atomic number. The elemental abundances are given relative to 
the abundance of silicon. 
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ordered by their Fe/O abundance ratio, deviate in a smooth monotonic fashion 
from the "mass unbiased baseline" element abundances 

Webber (1975) first noticed a systematic difference between SEP element 
abundances and the elemental abundances of the solar photosphere Webber 
showe*. that the element abundance ratio SEP/ photosphere was correlated with 
the element's first ionization potential Cook and Meyer have confirmed this 
behavior with their SEP abundance' averages Figure 5 2 2 from Cook (1981) 
shows the element abundance ratio SEP/photosphere plotted vs first ionization 
potential. The SEP(Cook) /photosphere ratio is of the order of 1 or larger, for 
elements with first ionization potential. 1, less than 10 eV (including Ni, Fe. Ms, 
Si) The SEP/photosphere ratio is less than one half for elements with 1 > 10eV 
(including C. N. 0. He. Ne. Ar) McGuire et al have noticed a similar behavior 
for their sample of flares which have low Fe/0 abundances Another possibility 
that has been suggested for the correlation between the SEP/photosphere 
abundance ratio and the first ionization potential, is a systematic error in the 
photospheric abundances, organized by 1 (see eg Mcguire et al 1979) 

Cook (1981) in analysis confirmed by Meyer (1981b) has shown that the 
SEP element abundance pattern is more similar to the element abundances 
observed for the solar corona and solar wind than to the photosphere Figure 
5.2 3, from Cook (1981), shows the ratios of the SEP(Cook) element abundanc es 
to the element abundances in the photosphere, corona and solar wind Cook 
has stated the hypothesis that the reason solar wind, coronal and SEP average 
abundances are similar, is that the solar wind and SEPs are accelerated coronal 
material and that coronal elemental abundances are significantly different from 
element abundances in the photosphere Further evidence supporting this 
hypothesis consists of the charge states of heavy Ions observed in the solar 
wind (Bame et al 1979 ), the charge slates c. low energy SEPs of « 10 to 800 
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KeV/nuc (Sciambi at at. 1977 and G'oekler at a l. 1976), 
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flgurt.' 5.2.2 - ( from Cook (1981)) SEP elemental abundances divided by 
solar photospheric or "Local Galactic" abundances of Meyer (1979) are 
plotted vs. the first ionization potential. 
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figure 5-2 3 - ( from Cook(196l)) The four flare SEP elemental abundances 
SEP(Cook) are compared to photospheric, coronal, and solar wind ele'— n- 
tal compositions. 
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and the charge stales indirectly inferred for SEPs of 0.5 to 5 MeV/nuc 
(O'GalL'ghe. ef al 1976) The most consistent explanation of the charge states 
observed is that matter accelerated to solar wind or solar energetic particle 
energies was once in collisional equilibrium with a gas that had a temperature 
of about 1-2 x 10* ' K. This temperature is about the same as observed in the 
corona (Athay 1976), and different from the temperature of the photospheric 
minimum, which is about 5 x 10 3 * K 

The pathlength of material that SEPs have traversed from acceleration to 
observation, as revealed by nuclear fragmentation, is consistent with the 
KUggest'on that the corona is the SEP particle source McGuire at al (1979) set 
a limit of 0 06 g / cm. 7 f. jm observations of the SEP *H/ H abundance ratio and 
Cook(l90l) reports upper limits tc B and F abundances that are consistent with 
the grammage reported in McGuire at al (1979) 

If the corona is the source of SEPs and the solar wind, what conclusions 
might we draw for Isotope observations? Since the difference between the pho- 
tosphere and ccrona or S> abundance averages seem to be organized by first 
ioiizatlon potential, the process that du.erenuates the elemental abundances 
of the corona from those of th' photosphere may be influenced primarily by 
whethf* a givr n irn is neutral or singly charged ir the environment where the 
differentiation takes place A way to explain the con elation of 
SEP/photosphere with i, is to assure* that elements with 1 > OeV are neutral 
and elements with 1 < "OeV a^t tingly charged, which aught imply that the 
efficiency of transport from photosphtra to corona for neutral atoms is 25% to 
40% of the transport efficiency foi ions (Meyer tb?-’b) 

Nl. Fe. kg. and Si all have I ^ 0eV and thus might '« singly charged in 
Meyer s picture of the transport process from photosphere to corona Ni. Pc. 
Mg. and SI. all have abo’*‘ the same ratio of SFP/photosphere elemental abun- 
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d&nces. If the mass of the Ion was Important to the transport process from 
photosphere to corona, then we might expect to see a different ratio or 
SEP/photcsphere for Ni and Fc than for Mg and Si. Cook(198l) observes about 
the same ratio of SEP/photosphere for Ni and Fe and for Mg and Si. Thus, a 
model that ties the differences in SEP/photosphere abundance to first ioniza- 
tion potential would have the isotopic composition of the corona and any SEP 
isotopic abundance average be the same as the isotopic composition of the pho- 

I 

toaphere. 

Table 5.2. 1 lists the SEP(Ccok) abundances and the abundance results 
obtained in section 4.3. Also Listed are abundance measurements by McGuire 
at al. (1979) in the energy ran go 6.7-15 MeV/nuc, for the time interval 
70 206:20 - 78 270:00, and a measurement by von Kosenvinge and Rearaes (1979) 
for day 70 267, in the energy range 6 7-12 4 MeV/nuc. Note that data from the 
flare onset is included in the von Rosenvinge and Reames, and McGuire at a L 
abundances. Cook (1981) has systematically excluded data in flare onsets from 
his SEP(Cook) averages Von Rosenvinge and Reames (1979) report a large 
abundance ra^.o variation for Fe/0 as a function of energy on day 78 267 They 
report Fe/O abundance ratios of 0 37 ± 0.02, 0.13 ± 0 01, and 0 067 ± 0 007 for 
energy ranges of 2 0-3 1 MeV/nuc, 3. 9-0. 7 MeV/nuc, and 6 .7-12 4 MeV/nuc, 
respectively. Figure 5.2.4 is a plot of the ratio of the abundances observed in 
the 78.266 flare to the SEP(Cook) average. The 78:286 flare element abun- 
dances show a relative depletion in the elements Ne.Si.S.Ca, and Fe, relative to 
SEP(Cook) The magnitude of the depletion is consistent with the flare to flare 
variations seen by Cook T ho disagreement between our measurement of Ne/0 
and He/0 and those of McGuire at al , and von Rosenvinge and Reames is 
greater than would be expected from a statistical fluctuation Some of the 
disagreement for He, Ne, and Fe may be due to the fact that McGuire at al , and 
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von Rosenvinge and Reames included the flare onset in their abundance meas- 
urements. We conclude that the abundance pattern of the 78:266 flare is con- 
sistent with that measured for other large flares. 

In figure 5.2.4 neon seems to bo more depleted than might be expected 
from the smooth flare to flare behavior of figure 5.2.3. To investigate whether 
neon was especially depleted in the 78:288 flare, for figure 5.2.5, we plotted 
Ne/O vs. Mg/O for the seven flares In Cook (lS8l) and two other flares in 
McGuire at a L (1979) not reported in Cook (1981). Relative to the other large 
SEP events, the 78 268 flare has a typical value for Ne/O and for Mg/O 
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SEP Dement Abundances 

(relative to the oxygen abundance) 


Table 5.2.1 - 



(1981) 

SEP(Cook) 

(1979) 
von Rose riving e. 

(1979) 

McGuire 


Dement 

Reames 

at a l. 

this work 

He 

72 ±3.3 

— 

40 ±3 

53 ± 2“ 

C 

0.47 ± 0.03 

0.40 ± 0.03 

0.49 ± 0 03 

0.47 ± 0 03 

N 

0.12 ± 0.01 

0.10 ± 0 01 

0.12 ± 0.02 

0.13 ± 0 01 

0 

= 1.0 

= 1.0 

= 1.0 

= 1.0 

Ne 

0.17 ± 0.01 

0.09 ± 0.01 

0.09 ± 0.01 

0.13 ±0.01 

Mg 

0 21 ± 0 02 

0.21 ± 0 02 

0.18 ± 0.02 

0.20 ± 0 02 

Si 

0.17 ± 0.01 

0.15 ± 0.02 

0.18 ± 0.02 

0.13 ±0.02 

Fe 

0.20 ± 0 01 

0.067 ± 0.007 

0.08 ± 0.01 

0.14 ± 0.03* 
0.07 ± 0 03 6 
0.03 ± D OS'* 
0.08 ± 0 02* 

Energy 

5-15 Mev/nuc 

8 7-12 4 MeV/nuc 

8.7-15 MeV/nuc 

fits 

5-50 MeV/nuc 

Time 

4 flare average 
19 77-1978 

78:287 

78:200:20 - 
78:270 00 

78:268:00 - 
78 272 00 


• 78:209-76 272 

• RANGE 0. 8-10.4 MeV/nuc 

e RANGE 1. 11 8-18 5 MeV/nuc 
4 RANGE 2, 18 5-24 3 MeV/nuc 

• Average RANGES 0-2 
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flgure 5.2.4 - A plot of the ratio of the elemental abundances observed in 
the 78:206 flare to the SEP(Cook) element abundance average. The open 
triangles are measurements by McGuire at al. (1979), the crosses are by 
von Rosenvinge and Reames (1979), and the open squares are for this work. 
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figure 5.2.5 - A plot of the element abundance ratio Ne/O vs. the ratio Mg/O 
for several flares. The open and filled circles are from Cook (1981). The 
flares used in the SEP(Cook) average are drawn as filled circles. The 
SEP(Cook) average is an open circle with a dot inside The open triangles 
are from McGuire at at (1979). The open square is from this work. 
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The flare to flare element abundca.ce variations are a relatively smooth 
monotonic function of charge, about the SEP average. The maximal variation of 
the Fe/O abundance ratio in the Cook(1981) flare sample is from 0.045 to 0.87, 
a factor of about 15. Suppose that Fe and O had charges characteristic of a 
temperature of 1.5x10® *K. Fe* 10 ® and 0*®. as calculated by Jordan (1909). 
Suppose that the enhancement of the Fe/O abundance ratio, {•(FW/0). was 
related to the charge to mass ratios of Fe and 0 by the relation, with a varying 
from flare to flare, 


log 




(5.2.1) 


q/m( 0) 

Then for {( Fk/ 0) =15. a = -3.9:. Then the enhancement of the “Ne/^Ne 
ratio, predicted by equation r ^ i: £( 2 *Ne/ 80 Nel = 1.45. A relation of this type 
predicts that the v Ne/ ab^darvve ratio vari: tions eye correlated with the 
Fe/O abundance ra io variations Thi „ign of the elect is such that ^Ne/^Ne 
would be enhanced if Fe/O was enhanced. In the 78 268 SEP event, the Fe/0 
abundance ratio is not enhanced, so we would not expect the ^Ne/^Ne abun- 
dance ratio to be enhanced by this mechanism 


Solar Energetic Particle Isotopic Composition 

The isotopic composition of SEPs has been measured directly by instru- 
mentation on spacecraft and indirectly by analysis of moon rocks Direct meas- 
urements of the isotopic composition of SEP events for elements heavier than 
helium, have also been reported by the Umversity of Chicago charged particle 
telescope on IMP-8 The Chicago group has achieved a r.m s nu.ss resolution of 
« 0.42 amu for neon events and « 0 56 amu for magnesium evems. The IMP-8 
instrument has a 750 n thick front detector, and thus has an energy threshold 
for two parameter analysis approximately equal to the HIST threshold for 


RANGE 4 events 
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To make an isotope abundance measurement of moderate statistical accu- 
racy, with such a high energy threshold, the Chicago group summed over a 
group of large solar flares. Dietrich and Simpson (1979) report 
“Ne/^Ne - 7.7 +2.3 —1.5, (equivalent to “Ne/^Ne = 0.13 ± 0.03), for a group 
of 7 flares, from July 3, 1974 to September 14, 1977. Dietrich and Simpson 
(1979) divide the seven flares into two groups. 3 are iron-rich, with Fe/0 = 
0.80-1.14, and 4 are iron-poor, with Fe/0 = 0.031 - 0.22. For the set of iron-poor 
flares, they measure ^Ne/^Ne = 0.11 ± 0 04, and for the set of Fe rich flares, 
^Ne/^Ne = 0.15 ± 0.06. They And no statistically significant disagreement 
between the “Ne/^Ne abundances that they report for the two groups of 
flares. The Dietrich and Simpson (1979) measurements would s‘ill be consistent 
with the predictions of equation 5.2.1. which would predict an enhancement of 
» 35% for the “Ne/^Ne abundance ratio in the Fe-rich flares, relative to the 
Fe-poor flares. Our measurement of Fe/0 in the 78:206 flare would classify the 
78:206 flare as an Fe-poor flare. 

There is a need for more measurements of SEP isotope abundances to 
determine whether the flare to flare isotope abundances also vary in a manner 
analogous to the SEP elemental abundance variations. 

Dietrich and Simpson (1981) report ^Mg/^Mg = 0.13 ± 0.04, for a group of 
10 flares, from July 3. 1974 to April 28, 1978 They do not report results for the 
zl Ne/ 20 Ne and ^Wg/^Mg abundance ratios, because those isotopes were not 
resolved in their data. 

Indirect measurements of SEP isotope abundance ratios have been 
reported, using moon rocks as collectors for the SEPs Black ( 1972a, b) has 
presented evidence for the existence of a neon component in lunar soils and 
meteorites that was implanted by solar energetic particles of 1-10 MeV/nuc. 
Black finds isotope abundance ratios for the component. neon-C, of 
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s He/ 4 He = 4.1 ± 1.0 x 10~\ “Ne/^Ne = 0.094 ± 0.003, and 30 Ar/ M Ar = 4.1 ± 0.8. 
Venkatesan et a l. (1901) studied etched feldspar separates of luDar soils and 
conclude that the long term average solar flare Ne composition is about 
“Ne/^Ne « 0.083. Yeniv and Marti (1981) find “Ne/^Ne = 0.002 ± 0.001 in the 
near surface layers of a lunar rock and conclude that *°Ne may be depleted in 
solar flares. Etique et at. (1981) And ‘^Ne/^Ne i 0.089 for solar flares, by etch- 
ing lunar plagioclase samples. All experimenters me as: ire a “Ne/^Ne abun- 

i 

dance that is lower than the solar wind values reported by Eberhardt et a l. 
(1970), but higher than the direct measurements in this work and by the Chi- 
cago group. 

The a He/ 4 He abundance ratio has been measured by a number of groups of 
experimenters. Anglin (1975) measures 3 He/ 4 He = 0.027 ± 0.005 in a 1969-1972 
flare s um. Hurford (1975) ha3 3 He/ 4 He = 0.009 ± 0.004 in an October 1972 to 
November 1973 flare sum. For the large SEP events analyzed by McGuire et a L. 
(1979), only upper limits, ranging from 0.008 to 0.06, for the s He/ 4 He abun- 
dance ratio were reported. For the large August 1972 series of events, Webber 
et at. (1975) report finite 3 He/ 4 He abundance ratios that vary with energy, 
ranging from about 0 001 at « 10 MeV/nuc, to about 0 015 at 60 MeV/nuc. 
They interpret their result as due to spallation in the solar atmosphere. 

For another class of SEP event, higher values of 3 He/ 4 He have been 
reported. *He-rich (Garrard et a l 1973) events are SEP events in which 
3 He/ 4 He is noticeably enhanced ( 3 He/ 4 He ^ 0.1) over average values measured 
in large SEP events, discussed in the ai>* *.e paragraph. The high abundance of 
3 He/ 4 He coupled with the low observed abundances for D/H and T/H, means 
that high energy nuclear spallation reactions cannot have caused the high 
3 He/ 4 He abundances (Garrard et at. 1973) Relative to the average SEP abun- 
dances measured in large events, 3 He-rich events have lower H fluxes, higher 
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4 He/H ratios and higher ratios of elements with charges Z> 6 to 4 He. *He-rich 
events have generally been observed as small SEP events (Hurford at a 1. 1975, 
Hovestadt el a l 1975, Mason at al. 1980). In two small *He-rich events, men- 
tioned in section 3.4, » 300 He events were observed, while for the 70 288 flare, 
over 10" He events were seen. 

3 He-rich events are of interest here because of the large isotope enhance- 
ments observed; Hurford et al. (1975) report 3 He/ 4 He ~ 8 for one event, or 
2 x 10 4 times the value reported for the solar wind. It is possible that a 3 He-rich 
event could occur simultaneously with a large SEP event, or be triggered by it. 
Despite the low flux levels of 3 He-rich events, the extremely large isotope 
enhancements might bias the isotope measurements in large SEP events, if a 
3 He-rich event occurred at the same time as a large SEP event. 

According to a model by Fisk (1978), selected ions can be preheated at the 
flare site, prior to the flare, when the ions satisfy am electrostatic ion cyclotron 
wave resonamce condition The conditions necessary for the Fisk model to 
apply are that i H a? 0 2. T,/ T t i 5, amd that p, the ratio of thermal energy 
to magnetic field energy, satisfies p^ 10 -3 , at the flare site. The resonance fre- 
quency lies between the H and He gyrofrequencies. amd only the 3 He gyrofre- 
quency satisfies the resonance condition with its first harmonic. 

Ions other than 3 He can also be preheated in this moael if the second har- 
monic of the ion gyrofrequency fit satisfies the equation (Mason et al 1980) 

i_ 

1.05 < < 1.19 + 0.13(m/*/nO 2 . (5 3 1) 

figure 5 3.1, modified from figure 10 in Mason et al (1980), is a plot of the fre- 
quency of the second harmonic of the ion gyrofrequency for the various charge 
states of abundant solar isotopes Selective isotope enhancements for heavy 
elements in the Fisk model may not be plausible First, the equation 5 3 1 reso- 
nance condition selects a wide band of frequencies, 
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figure 5.3.1 - (adapted from Mason at at. (i960)) A plot of one half of the 
second harmonic of the gyrofrequency for various charge states of several 
ions, as a function of 0*/ Oj* . 
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and different isotopes of the same element are separated by much less than the 
width of the allowed frequency band. Secondly, in a realistic model of the flare 
site, the charge state distribution for the heavy ions would be spread over 
several charge states due to inhomogeneities in temperature at the flare site. 
For neon the primary charge state varies from +6 to +10 as the temperature 
varies Trom « 5 x 10° *K to » 5 x I0 e • K, according to the charge state calcula- 
tion* of Jordan (1989). In particular, a large “Ne/®Ne enhancement could 
only be produced if most of the neon was in charge state +”. at a temperature 
of * 5 8 x 10° *K, and very little neon was in charge state +6. and if the upper 
limit to the enhanced gyrofrequencies was about 2 1^/0/* = 1.3. Our upper 
limit for the *He/ 4 He abundance ratio in the 78266 SEP event is 
3 He/ 4 He < 0.0028, measured in RANGES 2-5 

Thus we believe that the 78 266 flare is not described by this type of flare 
model. The narrow range of flare site temp matures required to produce a 
^e/^Ne enhancement and the absence of any measurable enhancement in 
the *He/ 4 He abundance ratio argue against the Fisk (1978) model for flare site 
preheating in this flare 

5.4. Solar Isotope Measurements 

Direct measu-ements of the solar isotopic composition are rare Some stu- 
dies have looked at the spectra of molecules in sunspots, which are a cooler 
environment than the photosphere Hall (1973) and Hall Noyes, and Ayres 
(1972) have investigated the infrared spectra of the molecule CO The abun- 
dance of ,3 C/ **C was found to be 0 94 ±01 of the "terrestrial" value. ,8 0/ ,8 0 
was measured as 1 08 t 0 35 of the terrestrial value With Cameron s (1981) 
values for the the above terrestrial abundance ratios, then 
,3 C/'*C = 0 0104 ± 0 0011 and 18 0./ ,8 0 = 0 0021 ± C 0007 
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Boyer. Henoux. and Sotirovski (1971) observed absorption lines due *.o the 
molecule MgH in sunspot spectra and obtained ra Mg/ 24 Mg = 0.125 t 0.030 and 
**Mg/ M Mg = 0.125 ± 0.038. Their results conflict with Branch (1970) who 
obtained “Mg^Mg^Mg of 60.20:20. Branch's work was criticized in lambert. 
Maliia, and Petford (1971), who set limits for the abundance ratios * 9 Mg/ Z4 M g 
and u Vg/ u iig. that they believe are consistent with the terrestrial abundance 
ratios of ^Mg^Mg^Mg of 80:10:10. 

Hall (1975) has measured solar 3 He/ 4 He abundance by fitting the shape of 
a helium spectral line feature at 10830 A, observed in emission in a pror^nence 
Hall obtains 4 ± 2 x 10 -4 for the *He/ 4 He abundance ratio. 

5l 5. Solar Wind Isotopic Composition 

The solar wind provides a sample of material from the sun. The isotopic 
composition of the solar wind has been measured by instrumentation flown on 
satellites, and by mass spectroscopic means, after collection on the surface of 
moon rocks, or on special foils deployed by the Apollo astronauts. For selected 
isotope ratios, the composition of the solar wind has been experimentally meas- 
ured to an accuracy of a few per cent. The major problem is to relate the isoto- 
pic composition of the solar wind to the isotopic composition of the sun 

The elemental and isotopic composition of the solar wind is not constant as 
a function of time. The He/H abundance ratio in the solar wind can vary on 
time scales from years to minutes, the magnitude of the He/H values reported 
have ranged from 8 1 x 10~* to 0.42 (Neugcbautr 19B1). The s He/ 4 He abun- 
dance ratio is also observed to vary O s iivie et a l (1980) find that 3 He/ 4 He is 
correlated with 4 He bulk velocity and varies from greater than 5 x 1CT 3 to less 
than 2 x 10“° Thus time variable mas.' fractionation effects occur in the solar 


wind. 
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The long term average abundances of elements and isotopes in the solar 
wind do not show such a large range of venation. Ogilvie and H^shberg (1974) 
And a solar cycle variation in the He/H ratio of only 0 01 t 0.01. from 0 035 to 
0.045 The abundance ratios, though relatively stable, may still be appreciably 
fractionated and different from the Sun's composition The sclar He/H abun- 
dance is not well known Heasly and Milkey (1978) measure He/ H = 0.1C ± 0 025 
in quiescent solar prominences. Ross and Alter (1976) have 
He/H = 0.063 ( + 0.037 - 0.023). based on Lambert's (1967) analysis of the pho- 
tospheric O/H abundance and the He/0 abundance from SEPs. 

If the He/H abundance in the outer parts of the sun has not been altered 
by nuclear processing, then it should be equal to the He/H ab'-mdance at the 
time that the sun formed Geiss and Reeves (1972.1981) set an upper limit of 
2 x 10"° for the D/H ratio in the material out of which the sun was formed This 
upper limit to D/H. according to a standard model of b.g-bang nucleosynthesis, 
(Wagoner 1973), sets a lower limit to the He/H abundance of He/H = 0 075, as 
the stellar destruction of He since the big bang is negligible 

A reliable value of He/H in the solar neighborhood is He/H = 0.101 t 0 005 
(Piembert 1975). obtained by analyzing H D regions of the Orion nebula. Thus, 
the sun probably has a value of He/H of about 0 1 and He/H is very likely to be 
greater than 0 075 The ratio of these estimates of the solar He/H value to the 
long term averages of solar wind He/ H ratios is at least 0 075/0 045 = 1 67. and 
probably greater than 2 Further evidence for helium depletion in the solar 
wind may come from the helium enhancements observed in the post-shock 
solar wind plasma, after the passage of a flare associated shock wave (Bame 
at ai. 1979. Hundhausen 1972). 

Models of the solar wind are not yet very successful at explaining the rea- 
sons for the fluctuating abundance ratios ’He/’He and He/H Geiss et al 
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(1970) have considered a model of the solar wind in which a proton fluid 
expands out from the solar corona, dragging the ions by dynamical friction. In 
that model ions are accelerated in the solar wind only if the proton fluid has a 
certain minimum flux, which is a function of the ion's charge and mass. For 
charge states characteristic of the solar corona, if 4 He is accelerated in the 
solar wind, then so are all other ions heavier than 'H. with the exception of D. 

Ogilvie at al. (1980) h»ve searched for a correlation, predicted by the 
model of Geiss «< al . , between the 3 He/ 4 He abundance in the solar wind and the 
H flux. They do not find such a correlation, but instead And a correlation 
between 3 He/ 4 He and 4 He flux. Schmidt at al. (1900) measure the speeds of H*. 
He* 4 , and 0 8 * in the solar wind and And v(He* + ) - u(H*) and v(0®*) - v(H + ) to 
be on average equal and in the range 0 to +60 km/ sec. Ions with velocities 
greater than hydrogen require a means of acceleration in addition to dynamical 
friction (Schmidt at al. suggest solar wind plasma waves) which could change 
Geiss et al. 's (1970) conditions for fractionation effects. 

We can conclude that a consistent theoretical and exper'mental framework 
to explain the helium element and isotope abunaance variations in the solar 
wind does not yet exist. In particular whether or not there is a difference 
between the long term solar wind 3 He/ 4 He abundance and that of the sun is not 
known. Geiss et al. (1972) measure 3 He/ 4 He = 4.3 ± 0.2 x 10~* in the solar wind 
and Hall (1975) obtains 3 He/ 4 He = 4 ± 2 x 10 -4 for the spectra of a solar prom- 
inence If the difference between solar wind and the sun's He/H abundances 
depended only on the charge to mass ratio, then the 3 Hc/ 4 He abundance might 
be « C;3% larger than on the sun. We based this estimate on equation 5.2.1 and 
an assumed depletion of the He/H abundance ratio in the solar wind relative to 
the Sun of a factor of 2. The data do not exclude this possibility. 


The solar wind noble gases have been extensively measured after collection 


on special foils and on moon rocks. Geiss ft oi (1972) averaged the results of 
several Apollo foil collection experiments and obtained 
3 He/ 4 He = 4 3 t 0 2 x 10”\ n Ne/ w Ne = 0.073 t 0 002, and 

^Ne/^Ne = 0 033 t 0 004. Kberhardt ft oi (1970,1972) analyzed trapped solar 
wind noble gases In lunar soils and obtain *He/ 4 H# = 3 7 t 0 1 x 10 4 and 
^Ne/^Nes 0 078 tOOOl for the present day solar wind There is a slight 
cLfference between the two sets of values, which is evidence that the lighter Iso- 
topes are slightly depleted in the moon rocks. 

Kberhardt ft ai (1972) report that for trapped Ar in lunar soils, the abun- 
dance ratio is ^Ar/^Ar = b 33 t 0 03, in agreement with terrestrial ^Ar/^Ar 
vnluos They report trapped Kr isotope abundances in lunar soils that are frac- 
tionated by about 0 hX per a mu relative to terrestrial values, but corsider 
solar wind and terrestrial Kr to be consistent within the accuracy of the 
method The lunar trapped Kr isotopic abundances also agree with the average 
value of trapped Kr In carbonaceous chondrites The agreement of Ar and Kr 
isotope abundances with terrestrial values argues against solar wind mass frac- 
tionation of the noble gases, unless all of the terrestrial Ar and Kr was 
deposited on the earth by the solar wind If there was solar wind isotope fractio- 
nation. then the terrestrial pattern would have to exactly match ♦ he solar wind 
pattern, to reproduce the experimental results 

The isotopes of Xe do not agree with the above picture (Podosek *978. 
Kberhardt nt oi 1972) Compared with terrestrial values, the trapped solar Xe 
isotope abundances show mass fractionation of about 3X per a nvu . difficult to 
explain at the same time as the Ar and Kr results The average isotopic compo- 
sition of carbonaceous chondrites agrees with solar trapped Xe for mass 
numbers 124-130 Xe mass numbers 132*138 for the chondntes's abundances 
show an excess compared to the solar trapped Xe. which has been attributed to 
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fission fragments However, the mass spectrum of the Xe is not that expected 
from common fission parents Xe then may question the internal consistency of 
the above picture that the solar fund has undergone minimal mass fractiona- 
tion 

Other elements have been identified on the surface of moon rocks as hav- 
ing possibly been implanted there by the solar wind Thero is evidence that the 
surface component of lunar nitrogen has had its '®N/ U N abundance ratio 
increase in the past by about 15-20% (Kerrtdge *t a i 1977, Becker and Clayton 
1975) For samples in which both ia N/ U N and w Ne/ w Ne were measured, the 
change in the l *N/ 14 N abundance is not correlated with the w Ne/ w Ne abun- 
dance Becker (1980) finds a weak correlation betwoon the l8 N/ u N abundance 

and ‘ a C/ !, C abundance, though the change in ia CY 1f C abundance is about ^ 

as largo os the change in the ,8 N/ M N abundance Kerridge oi (1977) find no 
convincing explanation of the variable lB N/ abundance ratio, and they con- 
sider both effects at the lunar surface and mechanisms that change the solar 
wind composition with time 

Epstein and Taylor (1972) have found that carbon in lunar soils has a total 
carbon content that is correlated with the hydrogen content Samples with 
high carbon content have ia C/ l *C enriched by about 5% over soils with a low 
carbon content They do not believe that they have positively identified the 
source of the isotopically heavy carbon, though cne of the possibilities that 
they mention is that the heavy carbon comes from the solar wind 

Thus, for carbon and nitrogen, incomplete physical understanding of the 
source of lunar carbon and nitrogen prevent a sure identification of the isoto- 
pic composition of solar wind carbon and nitrogen We can distinguish several 
possibilities For the first, we may suppose that n thor the ,0 N/ U N nor the 
IS C/ ■•C abundance ratios measured above accurately reflect the composition of 
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the solar wind. For the second, if either the carbon or the nitrogen isotopic 
abundance accurately reflects the solar wind, then that would be strong evi- 
dence against a simple linear mass fractionation process operating in the solar 
wind acceleration process. If the nitrogen results accurately measure tbe solar 
wind and solar wind compositional time variations, and are not a result of a 
changing sclar surface composition, then the problem becomes more complex 
If some unknown effect can cause the ,0 N/ U N ratio to vary and not oth er iso- 

I 

tope ratios, like ®Ne/ ro Ne and l3 C/ l *C, then we cannot predict the behavior of 
one isotope ratio by observing another isotope ratio The nitrogen results, like 
the Xe results, might tend to question the internal consistency of of a simple 
picture suggesting the absence of solar wind mass fractionation for ’sotope 
ratios in the solar wind. 

Zinner et a L. (1977) analyzed surface enhancements of Mg that may be due 
to the solar wind. They And isotopic abundances that are consistent with ter- 
restrial values within the experimental uncertainties estimated (Zinner 1 979) to 
be about 20%. 

bJ. Comparison of SEP and Solar Isotope Abundances 

The goal of our analysis is to compare solar and SEP isotope measure- 
ments That task is made difficult by the scarcity of direct solar isotope meas- 
urements. As discussed in section 5 4, the direct measurements of the ' sotope 
ratios 5 He/ 4 He, l3 C/ 12 C, ,fl 0/ lfl 0. ^Mg/^Mg, and ^Mg/^Mg. have uncertainties 
of 50%, 10%, 35%. 30%. and 30%, respectively We could not And direct solar 
values for the isotope ratios K Ne/ TO Ne and ls N/ l4 N Table 5 9 1 contains a list 
of measurements of the isotope ratios of the elements He, C, N. 0. Ne, end Mg. 
discussed in this chapter, along with our SEP measurements 

A source for indirect information about the isotopic composition of tne Sun 
is the composition of the carbonaceous chondnte class of meteorite As 
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discussed in Wasserburg at oi (1980), the isotopic composition of meteorites 
and terrestrial materials is cnaracterized by its homogeneity, the isotopic com- 
positions of any two samples of solar system material are likely to be the same 
to an accuracy of much better than a per cent. This was thought to be strong 
evidence for the view that the solar system at "birth" had a uniform isotopic 
composition, and that the small variations in isotopic composition seen today 
on the Earth are the result of various physical and chemical fractionation 
processes. Recently, as discussed in Wasserburg at at. (1980), evidence of pri- 
mordial inhomogeneities has been found, in the anomalous isotopic composition 
of parts of certain meteorites. Except for some coble gases, the maxim um 
difference between the anomalous and terrestrial isotopic abundances that has 
been seen is about 5-8%. 

The best accuracy that we report for our SEP isotope ratio results in tables 
4 2.1 and 5.6.1 is about 20%. For our purposes we can take as "solar” the ter- 
restrial and meteoritic abundances for the elements which are not noble gases, 
C, N, 0. and Mg We have used the terrestrial and meteoritic abundance compi- 
lation of Cameron (1981) as a source for these values, also listed in tables 4 2.1 


and 5.6.1. 
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Table 5.6.1 - 

abundance ratio value 

*He/ 4 He 4.3 ± 0.2 x lO" 4 

4.0 x 2.0 x 10~* 
3.7 t 0.1 x 10^* 

4.2 k 0 8 x 10^ 

1.3 iOBx 10^* 

0.027 t 0 005 
0 009 k 0.004 
•s 0 0026 

,3 C/ ,a C 0 0104 t 0011 

0.0111 

0 0095 + 0 0042 
- 0 0029 

,6 N/ ,4 N a + lOXto 

* - 5X of 
terrestrial 

0 0037 

0 008 + 0 010 
- 0.005 


source 

reference 

solar wind 

Geiss at at 
(1972) 

solar spectra 

Hall (1973) 

solar wind 
implanted in 
lunar soils 

Eberhardt at at 

(1972) 

meteorites 

(He-B) 

Mazor at at 

(1970) 

meteorites 

(He-A) 

Mazor at at 
(1970) 

SEPs 

Anglin (1975) 

SEPs 

Hurford (1975) 

SEPs 

this work 


solar spectra 

Hall (1973) 


Cameron (1981) 

SEPs 

this work 


solar wind 

Kerridge at at 

implanted in 

(1977) 

lunar soils 



Cameron ( 1 98 1 ) 

SEPs 

this work 
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Table 5.6.1 - 
(cooL) 


abundance ratio 

▼aloe 

source 

reference 

i»0/ >«o 

0 0021 * 0.0007 

solar spectra 

Han (1973) 


0.0020 


Cameron (1901) 


0.0015 + 0 0011 
- 0.0007 

SEPs 

this work 

•'Ne/^Ne 

« 

0.0028 ± 0.0001 

Ne-B 

Black( 1972a) 


0.0024 ± 0 0003 

solar wind 

Geiss e t al. 
(1972) 


0.0026 ± 0 0001 

solar wind 
implanted in 
lunar rocks 

Eberhardt et al. 

(1970) 


0 0030 


Cameron (1981) 


^ 0.014 

SEPs 

Mewaldt ct al. 
(1979) 


« 0.014 

SEPs 

this work 

“Ne/^Ne 

0.073 *0 001 

solar wind 

Geiss b t al. 
(1972) 


0.078 * 0.002 

solar wind 
implanted in 
lunar soils 

Eberhardt et al 
(1972) 


0 OQ0 * 0.001 

Ne-B 

Black (1972a) 


0.122 * 0.006 

Nc-A 

Pepin (1987) 

1 

1 

: 

0.122 


Cameron (1981) 

t 

t 

i 

i 

* 

* 

t 

i 

f 

0.13 ± 0 03 

0 1 1 * 0 04 
0 15 * 0 06 

SEPs 

SEIPs, Fe-poor 
SEPs. Fe-rich 

Dietrich and Simpson 
(1979) 

it 

t 

i 

i 

t 

0 13 + 0 04 
- 0 03 

SEPs 

Mewaldt at al 
(1979) 

♦ 

i 

t 

i 

L 

0.109 0 026 

- 0.019 

SEPs 

this work 
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Table 5.6.1- 
(codL) 




abundance ratio 

value 

source 

reference 


0.126 ± 0.038 

solar spectra 

Boyer at at. 

(1971) 


terrestrial ± 20% 

solar wind 
implanted in a 
lunar rock 

Zinner (1977) 


0 129 


Cameron (1981) 


0.15 + 0.05 
- 0 02 

SEPs 

Mewaldt et al. 

(1981) 


0.140 + 0.043 
- 0.025 

SEPs 

this worn 


2e Mg/ 84 Ms 0.125 1:0.030 

solar spectra 

Boyer et al. 
(1971) 

terrestrial t 20% 

solar wind 
implanted m a 
lunar rock 

Zinner (1977) 

0.142 


Cameron (1981) 

0.15 + 0.04 

SEPs 

Mewaldt et al 

- 0.03 


(1981) 

0.13 ± 0.04 

SEPs 

Dietrich and Simpson 
(1981) 

0.148 + 0 046 
- 0.028 

SEPs 

this work 


The noble gases are a more complicated story. For the elements He and Ne 
observed in meteorites, two primary components have been discovered, ‘'A”, or 
the "planetary" component and "B", or the "solar" component. He-A has 
s He/ 4 He « 1.5 x 10 -4 and He-B has 9 He/ 4 He » 3.9 x 10 -4 , while Ne-A has 
"Ne/^New 0.122 and Ne-B has ®Ne/ w Ne « 0.000 (Podosek 1970). "B", the 
"solar" component, is thought to have been implanted in the meteorites by the 
solar wind. This idea is strengthened by the analysis or trapped solar wind 

i 

gases In lunar soils, where ^Ne/^Ne = 0.070 ± 0.001 and 
3 He/ 4 He = 3.7 i 0.1 x 10 -4 (Eberhardt et a l. 1970) is obtained. It has been sug- 
gested that the difference between the present day solar wind neon measure- 
ments of Geiss et a l. (1972). with ^Ne/^Ne = 0.073 ± 0 002, and the two meas- 
urements of solar wind neon trapped in lunar soils and in meteorites may 
reflect mass fractionation that occurs in the trapping process. At argon the 
difference between the "solar" and "planetary" components becomes much 
smaller. In meteorites Mazor et at. (1970) have found that there are only small 
variations in 38 Ar/ 38 Ar 38 Ar/ 38 Ar varies from 5.3 to 5.5, as the 28 Ne/ 20 Ne com- 
position varies from 0.125 to 0.080. or from the "planetary" to "solar" com- 
positions. 

We want to know the isotopic composition of the Sun, for He and Ne; is it 
similar to the "solar" or "planetary” meteoritic compositions, or is it something 
else entirely? We will explore the views that the Sun is either "solar" or "plane- 
tary” in its noble gas isotopic composition, and examine the consequences of 
each 

The argument that the Sun has "solar” noble gas isotopic composition is 
supported by the ncble gas elemental abundance patterns in "solar" and 
"planetary” gases Relative to the "solar” gases, the "planetary” noble gas ele- 
mental composition pattern shows a depletion of about 4 orders of magnitude 
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for the light noble gases relative to the heavy noble gases (Podosek 1978). This 
elemental depletion could have been accompanied by isotope fractionation. 
The lighter isotopes of the "planetary" component of He and Ne are depleted 
relative tc the "solar" component. 

The close agreement between the isotopic composition of the "solar” 
meteoritic gases for the elements He and Ne. and the solar wind measurements 
of He and Ne, support the identification of the "solar" meteontic component as 
implanted by the solar wind. Then the Sun would have " solar"noble gas isoto- 
pic abundances if the solar wind was an unfractionated sample of the Sun. In 
section 5.3 we discussed the evidence for the absence of a simple linear mass 
fractionation bias in the long term average isotopic composition of the solar 
wind. In general the evidence was consistent with a lack of a simple linear mass 
fractionation bias, but the elements nitrogen and xenon were not explained by 
the simple picture. If there is no mass fractionation bias to the solar wind, then 
the solar isotopic composition for the elements helium and nson would be that 
measured by Geiss at ai. (1972), with s He/ 4 He = 4.25 ± 0 22 x 10~\ and 
^Ne/^e = 0.073 ± 0.002 Our SEP measurement of 

^Ne/^Ne = 0.109 + 0 026 - 0.019 disagrees with the solar wind measurement 
and we would then have to look for bias in the solar flare acceleration process 

The view that the Sun is made of "planetary" noble gases seems to require 
some type of mass fractionation in the solar wind acceleration process, for the 
elements helium and neon. That view assumes that the depletion of lighter 
noble gases in the "planetary" component of meteorites did not afTect the isoto- 
pic composition of the gases left behind in the meteorites The "planetary" iso- 
topic compositions for helium and neon then more closely reflect the isotopic 
composition of the Sun. than does the isotopic composition of the solar wind. 
This view is supported by Cameron (1981) who uses "planetary" isotopic compo- 
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sitions for his noble gas isotopic compositions. Mewaldt at a l (1979,1981), in 
an earlier presentation of this work, have argued that their SEP measurements 
support the "planetary" view. Their SEP measurement of the “Ne/^Ne ratio is 
much closer to neon-A than to neon-B. The other isotope abundance ratios 
that they measure for C, N, 0. and Mg. do not show any evidence for a simple 
SEP mass fractionation bias that would explain the "Ne/^e result 

The isotope abundance ratio results in this work do not represent a statist- 
ically significant difference between this work and Mewaldt et a l (1979, 1901). 
1 do not find their suggestion that the Sun is made of neon-A attractive, pri- 
marily because 1 find no convincing evidence of solar wind mass fractionation. 
We have seen in the previous section that Ar and Kr isotope abundances for the 
solar wind implanted in moon rocks agree with terrestrial abundance ratios 
This is a strong argument against a simple linear solar wind mass fractionation 
According to the theory of Geiss at al (1970), Ar requires a higher flux of H in 
the solar wind to be "puiled along" in the solar wind than does Ne, and so should 
be even more susceptible to mass fractionation. To me it is more plausible that 
an isotopic mass fractionation accompanied the observed depletion of the light 
noble gases in ihe format'on of the "planetary" noble gas component of meteor- 
ites 

In the following discussion we will compare our measurements with our 
best estimate of the Sun's isotopic composition We will take the position that 
the solar wind isotopic composition is identical to that of the Sun. The direct 
measurements of Geiss et al (1972) will be used for the He and Ne abundance 
ratios For C, N, 0. and Mg, we will use the meteoritic and terrestrial isotope 
abundances from the Cameron compilation (1901). 

Table 5 6 2 contains the comparison between our measurements (corrected 
for energy /nuc) and the solar values The "reduced difference SEP - solar" is 
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the difference between our SEP measurement and the solar abundance ratio 
given in the table, divided by the error in the SEP measurement. There is good 
agreement between our SEP measurements and the solar abundances, except 
for “Ne/*°Ne. The isotope abundance ratio for which the best solar wind and 
SEP measurements exist, shows a marked disagreement between the solar wind 
end SEP values. Our likelihood confidence interval calculations give a 2 5X pro- 
bability that the SEP isotope ratio could be as low as ze Ne/®Ne = 0 073, meas- 
ured in the solar wind 

We can investigate the possibility that a simple linear mass fractionation in 
the SEPs produced the distribution of isotope abundance ratios that we meas- 
ured. We define a normalized isotope abundance ratio for element Z as the 
ratio of our SEP isotope abundance ratio Tor element Z divided by the solar iso- 
tope abundance ratio for element Z. Figure 5.8.1 shows the normalized isotope 
abundance ratios for our measurements of ,3 C/ 1Z C, ' B 0/' 8 0, ^Ne/^Ne, 
*°Mg/ 24 Mg, and M Mg/ 24 Mg, plotted as a function of the ratio of the isotope 
masses. A,/ where Aj is the heavier isotope For a simple linear mass 
fractionation process, the enhancement, {{j .i). of isotope j relative to isotope i, 
is 


= a 




(5.6 1) 


If a linear mass fractionation process had operated on SEPs, the plot of 
reduced abundances vs Aj / ^ would be a straight line. The line s intercept 
would be a normalized abundance equal to 1 at A)/ A % = 1 The data for ele- 
ments C, O, and Vg seem to be consistent with an absence of SEP mass 
fractionation The Ne point does not follow the trend of the other data 
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Figure 5.6.1 - A plot of the normalized isotope ratio abundance 
(Ai /X i ')sn>/ T as a function of the ratio of the masses of iso- 

topes i and j. \/ Aj. 
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Mewaldt at at. (1979) investigated the possibility that the “Ne/^Ne abun- 
dance ratio was altered by nuclear spallation reactions, which might have 
occurred as the SEPs pass through the solar atmosphere. They found that 
nuclear spallation reactions in the solar atmosphere were much more likely to 
produce *'Ne than “Ne. Then the observed upper limit to the t, Ne/ 80 Ne abun- 
dance ratio of “'Ne/^Ne £ 0.014. implies that much less than 
0.014/0.109 = 13% of the “Ne events could have been produced by high energy 

I 

nuclear reactions as the SEPs propagated from the flare site to HIST. It is then 
not likely that nuclear spallation reactions by tne SEP particles have caused 
the difference between the solar wind and the SEP “Ne/^Ne measurements. 

Thus, the neon seems to be a puzzle If neon is fractionated in the solar 
wind, the fractionation process must be chemically selective, it must produce 
much more severe fractionation for neon chan for argon and krypton. Tne 
same consideration applies to SEPs If neon is fractionated in our SEP meas- 
urements, then it must be fractionated more than the isotopes of C, 0. and Mg 

There is a model of the solar flare preacceleration process that seeks to 
explain the values of “Ne/^Ne reported in solar flares m this work and by 
Dietrich and Simpson (1975) as a solar flame preacceleration enhancement of a 
solar wind - like neon composition for the Sun. 

Mullan and Levine (1581) have introduced a theoretical model that seeks to 
explein the abundance enhancements observed in solar flares They model a 
flare preacceleration process that takes place inside a collapsing neutral 
current sheet, in a mag netic flux tube loop They assume that in such a struc- 
ture the maximum temperature of the loop may regulate itself 30 that the col- 
lapse time scale, f^. in only slightly smaller than the proton Coulomb less tirr 
tp. In such a situation, ions with atomic charges, are prv heated along 

with the protons bee anse they gain energy from the collapsing magnetic field at 
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a rate, l / tg, and lose energy at a rate, tj. auch that tj < t p < t g . The popula- 
tions of preheated ions might be converted to element abundance enhance- 
ments in SEPs by a SEP acceleration process that accelerated only particles 
above a certain energy threshold, for example. 

The Mullan and Levine model predicts that the SEP abundances for charges 
in the regions around Si and Fe will be enhanced compared to elements in the 
charge regions S-Cl and C-0. The abundance enhancements in this model may 
vary depending on assumptions made about conditions at the flare site. Mullan 
and Levine do not propose specific models for the possible flare to flare abun- 
dance variations. 

Mullan (1982) has applied Mullan and Levine's (1981) model to SEP isotope 
enhancements. Mullan finds enhancements of neutron rich isotopes of various 
elements for certain values of e in the relation 

( 5 . 02 ) 

that defines the relation between the magnetic collapse time and the Coulomb 
energy loss time. Mullan introduces a measure of inhomogeneity into trie 
preflare conditions at the flare site by integrating the fractions of various ions 
that are preheated over a distribution of t values. For a distribution for r. that 
is constant for |e| < 6c and zero for |e| ^ <5e. and for 6c » O.CJ, Mullan finds 
that the preheat mechanism would enhance the ”Ne/ Z0 Ne ratio more than for 
other measured isotope ratios 

The theory predicts that enhancements in the abundance ratio ,3 C/ l2 C will 
be produced for the same flare site conditions that produee the az Ne/ za Ne 
enhancements. Future high precision measurements or correlations between 
1S C’/ *®C end M Ne/®°Ne would be strong evidence for the Mullan and Irvine, and 
Midlan models Unfortunately, the Mullan and Levine theory has too many free 
parameters for a precise element by element comparison with SEP measure- 
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meats. 

There is a need Tor more measurements. SEP measurements of isotope 
abundances should be made for more flares, to assess the extent of variability 
of the isotope abundances. Difficult measurements of the present day soLar 
wine isotopic composition for C, N, 0, Mg, and other common elements that are 
not noble gases, should be attempted. 

5.7. Spectral Correction to Isotope Abundance Ration 

A subsidiary issue that has beeu leTt unmentioned so far is the question of 
making corrections to the isotope abundance r*tio^ based on considering the 
spectra as functions of energy /nuc or total momentum There arc two possible 
ways to discriminate between the two hypotheses One method to discriminate 
between the energy /nuc and total momentum isotope abundance corrections, 
is to plot isotope abundance ratios vs. energy and correct in turn for 
ene^gy/nuc and total momentum If one hypothesis produced an isotope abun- 
dance ratio that w»* approximately constant as a function of energy, and the 
other hypothesis did not, then a choice could be made Unfortunately, marginal 
statistics would make such a test inconclusive when applied to the isotope 
abundance ratios measured for the 78 266 flare. 

A second method considers the L.ntope abundance ratios integrated over 
all energies The ratio of the energy /nuc correction to the total momentum 
correction, for the 26 Mg/ 84 Mg abundance ratio in RANGES 2-4, is 0 858/1 444 -= 
0.53. a considerable difference Table 5.7.1 lists the isctope abundance ratio 
resuM for the 286 flare corrected for the totcu momentum functional depen- 
dence. Also listed in table 6 7 1 are the ratios of the correction factors for total 
rr^rnenturn to the correction factors for energy/nuc Table 5 7 1 lists the fac- 
tors by which to multiply the isotope ratio results of section 4 2 in order to 
obtam isotope abundance ratios corrected for total momentum or equivalently, 


- 170 - 

ngidtty. 

The reduced differences of SEP minus solar isotope abundance are all posi- 
tive, when we correct for total momentum and much larger than would be 
expected from random statistical fluctuations Then we conclude that the 
spectra of the heavier Isotopes are better represented by a spectral shape that 
; is a function of energy /nuc. 

&lB. Conclusions 

We have measured Isotope abundance in a large solar flare for the ele- 
ments of He. C. N, 0. Ne, and Mg The element abundances In the flare are not 
very different from the SEP average of Cook (1981). The 78 280 flare is slightly 
depleted In elements heavier than oxygen and in iron, relative to the four flare 
average of Cook(198l) The SEP isotope abundances observed in the Hare agree 

with terrestrial and solar isotope abundances, with the exception of neon The 

I 

neon isotopic abundance we measure for the flare. 
**Ne/*°Ne = 0 109 + 0 026-0 019. is different from the neon isotope abundance 
measured by Geiss ef ai (1972) for the solar wind. OT Ne/ w Ne = 0 073 t 0 002 
The neon isotopic abundance for this Hare is the same as ’’planetary neon”, or 
neon- A. which Cameron (1981) has adopted as a solar system value 

Separate arguments appear to rule out simple linear mass fractionation in 
the solar wind and in our SEP measurements as the cause of the discrepancy in 
the composition of the apparent compositions of these two sources of solar 
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SKP and Solar Isotope Abundances 
(SEP corrected for energy /nuc) 


Table 5.6.2 - 


reduced* 


Ratio 

solar 

abundance 

source 

SEP 

measurement 

difference 
SEP - solar 

*H&/ 4 He 

4.3 x 1(T* 

Geiss at a l 
(1972) 

< 0 0020 

— 

1S C/ 18c 

00111 

Cameron 

(1980) 

0 0095 + 0 042 
- 0.029 

•04 

ion/ U N 

00037 

H 

0 008 + 0 010 
- 0 005 

+0 9 

1Bq/ i*o 

0.0020 

II 

0 0015 + oorn 
- 0.U»'07 

-0.5 

“Ne/“he 

0.073 

Geiss at at 
(1972) 

0 109 + n C26 
- 0.019 

+ 1.9 


0.129 

Cameron 

(19P0) 

0 148 + 0 046 
- 0 028 

+0 8 


0.142 


0 148 + 0 043 
- 0 025 

+0.3 


* The reduced difference between the SEP and solar 
abundances is the difference between the STEP and solar 
abundances divided by the error bar for the SEP 
measurement 
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Isotope Abundances Corrected for Total Momentum 
Table 5 . 7.1 - 


Ratio 

(total momentum 
correction)/ 
(energy /nuc correction) 

SEP 

abundance ratio 
corrected for 
momentum 

reduced* 
difference 
SEP - solar 

ia c/ i*c 

1.390 

0.0132 + 0.0058 
- 0.0040 

+0.5 

, 8 N/ 14 n 

1.356 

4 

0.010 + 0.014 
- 0.005 

+ 1.2 

ia 0/ 18 0 

1.850 

0.0025 + 0.0019 
- 0.0012 

+0.4 

“Ne/^Ne 

1.611 

0.175 + 0.043 
-0.030 

+3.4 

“Mg / 24 Mg 

1.302 

0.192 + 0.060 
- 0.034 

+ 1.9 

“Mg / 24 Mg. 

1.607 

0.249 + 0.073 
-0.042 

+2.0 


• The reduced difference between the SEP and solar 
abundances is the difference between the SEP and solar 
abundances divided by the error bar for the SEP 
measurement. 


,1 
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Appendix A - Hodoscope Efficiency 

To calculate spectra we must be able to compute the efficiency, or the per- 
centage of a given type of particle eli mina ted, for various restrictions placed on 
the data set. In detectors Ml and M2, the "signal defect" problem, discussed in 
section 3.5.1, makes it difficult to determine the efficiency of hodoscope coin- 
cidence requirements directly from flight data. We identify elements in the low 
resolution data set by two parameter analysis. When we must use the matrix 
detector energies in the two parameter analysis, for RANGES 0 or 1, the " signal 
defect" problem could cause us to mlsidentify elements. For example, particles 
near the carbon track may be carbon nuclei or they may be oxygen nuclei, with 
a poise height defect. To be sure that we identify the elements correctly, we 
must require that only one matrix detector strip be triggered in each plane of 
the hodoscope Then to calculate spectra for each element anu RANGE, we need 
to know what fraction of the total events have single strip matrix detector 
triggers. 

If the percentage of events with double strip triggers varies with energy 
deposited in the matrix detectors, then analysis of flight data from RANGES 0 
and 1 to determine those percentages would involve a complicated deconvolu- 
tion procedure with simultaneous fits to element abundarces, energy spectra, 
and double strip fractions as a function of energy. Flight data from RANGES 2-8 
can be used over a portion of the necessary energy range, but suffers from sta- 
tistical limitations. Hius we chose to compute the double strip efficiency f~om 
cadibration data taken at the Bevalac in Berkeley. The same silicon detectors 
were used for Mi and M2 in flight and for the Bevalac run We looked for a rela- 
tion between the probability of a double strip event in either plane of a matrix 
detector and the energy deposited in that detector. 

An “Fe beam at « 600 AfeF/nuc was incident on a target of polyethylene 


i 
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(n(CHa)) The 8 ®Fe beam and the products of nuclear reactions in the target 
stopped in HIST. The same procedure was followed with a ,0 N beam at 
^ 600 HeV/nuc . As a result we collected data on nuclei with charges Z = 7 to 26 
that primarily stopped in the back detectors (D5 to D8) of HIST. 

We investigated the hypothesis that the fraction of events with double strip 
triggers in Ml or M2 depended on the total amount of energy deposited in Ml or 
M2. Three sets of * alibration data vfere considered, ^e ions that stopped in 
RANGES 2-0, “Si. “Si and ^Si ions that stopped in RANGES 4-0, and ,8 N ions that 
stopped in RANGES 4-0. Each RANGE was divided up into a number of energy 
intervals and a histogram or the Ml energy was accumulated for the single strip 
events m each interval 

Figure 3.8 la is an example of a histogram of the Ml energies for events 
that stopped in one of the RANGE 4 energy intervals, for events with single strip 
triggers. Figure 3 0.1b is a histogram of events that stopped in the same RANGE 
4 energy interval, and bad double stnp triggers in the "v plane” of Ml. The 
mean energy deposited in Ml was computed from the histogram for the single 
strip events The energy deposited in M2 is not read out of HIST, for RANGES 4- 
8 We calculated the energy deposited in M2 corresponding to the mean Ml 
energy, using a range-energy relation with the Barkas and Berger correction of 
section 2 6 4 

For each energy the percentage of double strip triggers in either plane of a 
matrix detector was computed Figure 3 8 2 shows the results plotted vs 
energy deposited in M2 There is a smooth trend with energy and the values 
derived using different charges agree in regions where the data sets overlap 
Data poiuts from oxygen RANGE 2 and 3 events observed in flight are also plot- 
ted. and agree with the calibration data The solid line in figures 3 8 2 is a fit to 
the calibration data, consisting of three line segments specified by six free 
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parameters. A quantitatively similar graph and fit were obtained for detector 
Ml. 

The functional form of the fit to the calibration data is given by the func- 
tion Photo { detector . £**), where "detector” specifies Ml or M2, and is the 
energy in MeV deposited in the specified detector. The highest energy for which 
this correction is used would correspond to a z8 Si nucleus that just stops in 
50 fi of silicon, the nominal thickness of Ml or M2, and is about 140 Me V. 

The matrix detectors are constructed with » 3/4 mm. wide metal charge 
collection electrodes or strips and a»l/4 mm. gap between the strips When a 
particle enters the gap between two strips, the charge divides and part goes to 
each strip. We might expect that as the energy deposited in the matrix detec- 
tor gap increases, it would be more likely for the amount of charge in both 
strips surrounding the gap to be above the threshold for triggering. We then 
expect the double strip probability to rise s'^wly with the amount of energy 
deposited in the matrix detector. The maximum probability for a doable strip 
in either plane of a matrix detector with a 1/4 mm gap would then be 
«(1 - 0.75) 2 * 44%. 

We do not understand the reason fcr the fall off in double strip probability 
above 100 MeV, but it is observed to be quantitatively similar for detectors Ml 


and M2. 
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figure 3.8. 1 a- b - Histograms of the energy measured Ln detector Ml, for 
“Fe events that stop in a restricted energy range of RANGE 4. Figure 
3.8.1a is for single strip triggers and figure 3.8.1b is for double strip 
triggers. 
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figure 3.8.2 - The probability of a double strip event in detector M2 is plot- 
ted as a function of the energy deposited in detector M2. Crosses are ,8 N 
calibration data in RANGES 4-8. solid squares are “Si. **Sl, and “SI calibra- 
tion data from RANGES 4-8, and the solid triangles are “Fe calibration data 
from RANGES 2-8. The open circles are oxygen flight data from RANGES 2 


and 3. 
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